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' We investigate the convergence properties of optimized perturbation tiieory, or linear S expansion 

(LDE), within the context of finite temperature phase transitions. Our results prove the reliability of 
W), these methods, recently employed in the determination of the critical temperature Tc for a system 

^ ' of weakly interacting homogeneous dilute Bose gas. We carry out the explicit LDE optimized 

' calculations and also the infrared analysis of the relevant quantities involved in the determination of 

Tc in the large- A'^ limit, when the relevant effective static action describing the system is extended to 
0{N) symmetry. Then, using an efficient resummation method, we show how the LDE can exactly 
reproduce the known large- A'^ result for Tc already at the first non-trivial order. Next, we consider 
the finite N = 2 case where, using similar resummation techniques, we improve the analytical 
^ ■ results for the nonperturbative terms involved in the expression for the critical temperature allowing 

ty^j ' comparison with recent Monte Carlo estimates of them. To illustrate the method we have considered 

I ! I a simple geometric series showing how the procedure as a whole works consistently in a general case. 
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O ; I. INTRODUCTION 

, Scalar field theories are extremely important in the study of symmetry breaking and restoration in different branches 
^ ' of physics such as cosmology, particle physics and condensed matter physics, where they may represent inflatons, Higgs 
\ particles, quark-anti-quark bound states. Cooper pairs, bosonic atoms and molecules, respectively. In most cases the 
' vacuum expectation value of those scalar fields represents an order parameter that signals phase transitions associated 
to symmetry breaking/restoration P]. 

In general, one important problem we have to deal with when studying phase transitions in field theory regards 
^vq the reliability of perturbation theory and its eventual breakdown. In this case, a nontrivial problem arises since 
C3 I nonperturbative methods must be used. This is the case in those physical situations involving a second order, 
' or weakly first order phase transition, where we have to consider the problem of infrared (IR) divergences that 
^2 i become progressively more important as one approaches the critical temperature, from above or below, and that will 
H ■ unavoidably spoil any perturbative attempt to compute relevant quantities there. In those situations we must find 
appropriate methods to take into account the large IR corrections, present in the form of large field fluctuations. 
There is a variety of nonperturbative methods that can be used in order to account for these corrections, including 
the recent dynamical Boltzmann-like approach that deals directly with the large field fluctuations 0|- At the same 
^ ; time, the most common methods, at equilibrium, try to resum the leading IR corrections. This happens for example, 
ILJ ' through the use of e-expansion techniques in order to compute corrections to the critical exponents that control the 
. !^ I singular behavior of physical quantities near the critical point (as it is familiar from the theory of critical phenomena 
^ ' [3)i the large- method and other approaches (for a review, see Ref. 0). 

?H I An issue that has attracted considerable attention recently and that is associated to the perturbation theory 
. . r breakdown problem is the study of how interactions alter the critical temperature (Tc) of Bose-Einstein condensation 
(BEC). Due to its nonperturbative nature, this is clearly a non-trivial problem. On the other hand, studies related 
to the BEC Tc problem are particularly important nowadays due the recent experimental realization of BEC on 
dilute atomic gases (for reviews, see for instance, Ref. 5]). The experimental achievement of BEC has induced many 
theoretical investigations which make use of methods developed to treat finite temperature quantum field theories. 
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At the same time, due to the high experimental precision with which the parameters may be tuned, BEC experiments 
provide an important laboratory to test many methods as well as models developed to treat those theories (see, e.g.^ 
Ref. 6]). 

The studies concerning the equilibrium properties of BEC can be addressed by means of a non-relativistic effective 
theory described by a complex scalar field. In the dilute limit, which is the regime involved in those experiments, only 
two body interactions arc important ^5] and one may then consider the following U{\) invariant finite temperature 
Euclidean action 



Se = dr y'd3^|v'*(x,T) (^^-^V^) V(x,T)-MV'*(x,r)V(x,T) + ^[V;*(x,T)^(x,r)]^| , (1.1) 

where, in natural units, f3 is the inverse of the temperature, fi is the chemical potential and m represents the mass of 
the atoms. At the relevant low temperatures involved in BEC the internal degrees of freedom are unimportant and 
this can be taken as an effective model of hard core spheres with local interactions for which a represents the s-wave 
scattering length. 

The field ip can be decomposed into imaginary-time frequency modes 'il'j{x,ujj), with discrete bosonic Matsubara 
frequencies Wj = 2nj/P, where j is an integer. Near the transition the chemical potential becomes very small as 
compared to the temperature {\fi\ << T) and, since the correlation length tends to infinity, the distances are large 
compared to the thermal wavelength A = yjiirfi/m. Therefore, the non-zero Matsubara modes decouple and one is 
left with an effective action for the field zero modes (j = 0) given by {Jl 



where -00 stands for the field's zero mode. Recently, Arnold, Moore and Tomasik have argued that when naively 
going from the original action {Se) to the reduced action (S'sd) by ignoring the effects of non-zero frequency modes 
one misses the effects that short-distances and/or high-frequency modes have on long-distance physics. For the critical 
temperature of condensation as a function of the density (Tc(n)), at second order, these effects can be absorbed into 
a modification of the strengths of the relevant interactions which means that one should consider the more general 
form for the reduced effective action Eq. 1)1. 2|l 



5eff[0O,0o] =P j d^X jVo ("^^^V" - ^0 + [0o*^o]' + O [^^VolV^P, (V-^V")'] | + /J^Vacuum , (1-3) 

where is the wave-function renormalization function, fi^ incorporates the mass renormalization function, Za 
incorporates the vertex renormalization function and i^vacuum represents the vacuum energy contributions coming 
from the integration over the nonstatic Matsubara modes. The O [-i/'gV'ol V'f/'oP, (V'oV'o)"^] terms represent higher order 
interactions in the zero modes of the fields. As shown in Ref. these terms will give contributions to the density of 
order-a"^ and higher and therefore do not enter the order-a^ calculations considered here. By matching perturbative 
order-a^ results obtained with the original action Se and the general effective action ^off, the authors of Ref. @ were 
able to show that the transition temperature for a dilute, homogeneous, three dimensional Bose gas can be expressed 
at next to leading order as 



,ln 



(1.4) 



where Tq is the ideal gas condensation temperature, Tq = I-k jm [7i/C(3/2)]^^'^, n is the density, C,{x) is the Riemann 
Zeta function and ci, d-^ and are numerical coefficients. A similar structure is also discussed in Ref. j^. As far the 
numerical coefficients are concerned, the exact value, — — 647rC(l/2)C(3/2)~^/^/3 ~ 19.7518, was obtained using 
perturbation theory Q. On the other hand, the other two coefficients, c\ and Cj, are sensitive to the infrared sector 
of the theory and consequently cannot be obtained perturbatively, but they can, through the matching calculation, 
be expressed in terms of the two nonperturbative quantities k and 7?. which are, respectively, related to the number 
density (i/'oV'o) and to the critical chemical potential /ic, as shown below. The actual relation between the two 
nonperturbative coefficients and these physical quantities is given by 



ci = -1287r3[C(3/2)]~4/3K ^ 



(1.5) 
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and 



647r 



C(l/2)[C(3/2)]-5/3inC(3/2), 



(1.6) 



where i^i Eq. is 



= 327r 



ln(1287r3) 



C(l/2)- 



^ [C(l/2)]^ 



(1.7) 



with i4r2 = —0.13508335373. The quantities n and TZ are related to the zero Matsubara modes only. Therefore, they 
can be nonperturbatively computed directly from the reduced action ^eff which, as discussed in the numerous previous 
applications, can be written as 



Srb — 



(1.8) 



where <j) = {(pi,(p2) is related to the original real components of ipo by ipo{x) = yJmT jZ^ [0i(x) + i(?!)2(x)], rbaro = 
ImyuzjZ^ and u = ASnamTZa/ Z^'^ . The vacuum contribution appearing in Eq. (|1.3|l will not enter in the specific 
calculation we do here and one has been omitted it from Eq. p.8|l . In the large- TV limit considered in the first part 
of this work, and also in Refs. the field (j) in Eq. H1.8|) is formally considered as having N components (0i, 

i — 1,. . . ,N). In this case, the Bose- Einstein condensate effective action Eq. (|1.8|) is the N — 2 special case of the 
general 0{N) invariant action. 

The three dimensional effective theory described by Eq. (|1.8|l is super-renormalizable and requires only a mass 
counterterm to eliminate any ultraviolet divergence. In terms of Eq. H1.8|l . the quantities k and TZ appearing in Eqs. 
(frB|l - (fTTjl are defined by 



(1.9) 



and 



(1.10) 



where the subscripts u and in Eq. H1.9|l mean that the density is to be evaluated in the presence and in the absence 
of interactions, respectively, and S(0) is the self-energy with zero external momentum. Since these physical quantities 
are dependent on the zero modes their evaluation is valid, at the critical point, only when done in a nonperturbative 
fashion. As discussed in the next section, the relation between Tc and S(0) comes from the Hugenholtz-Pines (HP) 
theorem at the critical point. 

Eq. 11.4|l is a general order-a^ result with coefficients that, therefore, depend on nonperturbative physics via k and 
TZ. In principle, to evaluate these two quantities one may start from the effective three-dimensional theory, given by 
Eq. (|1.8|l . and then employ any nonperturbative analytical or numerical technique. 

When quantum corrections are taken into account, the full propagator for the effective three dimensional theory 
reads 



G{p) = b2+r + S,en(p)]- 



(1.11) 



where represents the three momentum and Ercn(p) represents the renormalized self-energies. At the transition 
point (p^ — 0), the system must have infinite correlation length and one then has 



[G(0)]-i - [r, + S,e„(0)] = 



(1.12) 



This requirement leads to the Hugenholtz-Pines theorem result Vc — — Srcn(O). Since Tc is at least of order u it 
would be treated as a vertex in a standard perturbation type of calculation in which G{p) = 1/p^ represents the bare 
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propagator. This shows that perturbation theory is clearly inadequate to treat the BEC problem at the transition 
due to the presence of infrared divergences. One must then recur to nonperturbative methods like the numerical 
lattice Monte Carlo simulations, the analytical or the linear S expansion (LDE) adopted in this work (see 
for instance Ref. [l^ for earlier works on the method). The problem is highly nontrivial since the Hugenholtz- 
Pines theorem automatically washes out all momentum independent contributions, such as the one-loop tadpole 
diagrams, which constitute the leading order of most approximations. In practice, this means that the first nontrivial 
contributions start with two-loop momentum dependent self-energy terms. However, having reduced the original 
model, Eq. to the effective three-dimensional one, Eq. H1.8|l . makes it easier to tackle those contributions since 

one no longer has the problem of summing over the Matsubara's frequencies, which is a hard task when the number 
of loops increases. 

Recent numerical Monte Carlo applications 0, 0| have predicted values for ci which are close to 1.30 ^. On 
the other hand, some analytical applications have predicted values such as ^ 2.90 obtained with a self-consistent 
resummation method ~ 2.33 obtained with the 1/N expansion to leading order 0| and ~ 1.71 obtained with 
the same expansion to the next to leading order The LDE was first applied to order-J^ producing ci ~ 3.06 

[TtI I. Recently, the calculation has been extended to order-(5'' with the results ^ 2.45 at order-(5^ and ~ 1.51 at 
order-(5* The coefficient C2 was evaluated with Monte Carlo techniques and the predicted value obtained 

from those simulations is 75.7 ± 0.4. This quantity was also analytically evaluated with the LDE in Ref. where 
the encountered numerical values are ~ 101.4, ~ 98.2 and ~ 82.9 at second, third and fourth orders respectively. An 
order-(5^ application to ultra-relativistic gases has also been performed 19J. The LDE has been especially successful 
in treating scalar field theories at finite temperatures |20l |21| as well as finite temperature and densitv [221 . Several 
different applications performed with the LDE are listed in Ref. Recently, Braaten and Radescu j23j| have also 

used the LDE, with different optimization prescriptions, to evaluate Tc both at large and finite N limits, while Kleinert 
[23| has used the variational perturbation theory, which is a variation of the LDE, obtaining the value ci ~ 0.91 ±0.05. 
Like we do here for the finite N case, he has considered up to order-J"* contributions which include five loop diagrams. 
However, none of those authors considers resummation techniques to accelerate convergence and do not evaluate the 
coefficient C2, that it is also computed in the present work. 

As far as the application of the LDE to the determination of the BEC transition temperature is concerned, since 
the first papers making use of the LDE method to this problem Il8l |. an important question was raised and 
remained unanswered. This question regards the convergence properties of the method in this application for the 
BEC problem, which is in fact related to the convergence of the method in critical theories in general. Actually, this 
is a timely and important question regarding the applications of the LDE in field theories, since the first efforts have 
been concentrated mostly in the anharmonic oscillator problem at zero temperature, where rigorous LDE convergence 
proofs have been produced |2^ |2^ |23, ■ The extension to the finite temperature domain was also considered by 
Duncan and Jones p9l | , who used the anharmonic oscillator partition function. Only very recently a convergence proof 
has been extended, but for a particular perturbative series case, to asymptotically free renormalizable quantum field 
theories at zero temperature [s^l- Here, as in a previous paper from us, Ref. j^, our interest is to probe convergence 
in the vicinity of a phase transition, such as for the Bose-Einstein condensation problem presented above. This study 
should also settle the questions regarding the correctness of our original LDE applications (17i . ilS] , since by showing 
convergence in the large- iV case we also establish the reliability of our finite N results, as originally studied in those 
investigations. The LDE convergence in the large- iV extension of the BEC problem has also been recently addressed 
by Braaten and Radescu ■ The differences between their approach and ours will be discussed in some detail in the 
text. 

The literature shows [28L l3^ l3^ that in most applications it is possible to establish simple relations between the 
LDE and other nonperturbative methods already at order-J where one-loop diagrams are present. In fact, one can 
show that in those cases the LDE either exactly reproduces 1/A'^ results or produce very close numerical results. Here, 
the BEC problem poses an additional difficulty since, as discussed above, the first non-trivial contributions start at 
the two-loop level in the self-energies. As we shall see, it is not easy in this case to establish simple analytical relations 
for the quantities being computed, like those given e.g. in Refs. psL 133. l33l |. and the problem must be treated 
differently. However, as we are going to show in the coming sections, our numerical results, improved with an efficient 
resummation technique exactly converge in the large- limit and seems to also converge in the arbitrary N case. 

This paper is organized as follows. In Sec. II we briefly recall the LDE method and present the interpolated 
version of the action Eq. H1.8|l to be studied throughout the paper, following the recent applications performed in 
Refs. In Sec. HI, we carry out the formal evaluation of in three different ways. The first is the usual 

order by order type of calculation, which is familiar from perturbative calculations, and is in fact the only possible 



^ See Ref. II6II for an extension of these works to the 0(1) and 0(4) cases. 
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one for the realistic finite N case. The second uses the type of resummation famihar from nonperturbative methods 
such as Hartree and the 1/N approximations. These two procedures generate two series for the large- limit of (0^)u 
(in which case (j) in Eq. H1.8(l is extended to N components) whose coefficients, which are numerically obtained, can 
be usefully compared as a cross check, helping to establish the numerical reliability of the finite N series. In the 
same section, we also consider the asymptotic infrared and ultraviolet behavior of the series which, as we shall see, 
is a very useful approximation allowing at the same time a fully analytical analysis. In Sec. IV we first examine 
the LDE optimized perturbation procedure in a general case, considering a simple geometric series in order to get 
insight regarding the convergence structure of the optimal results. Then, the large- BEC series are optimized and a 
resummation technique which accelerates convergence is introduced. By taking the large- TV result for ci (obtained in 
Ref. a reference value we proceed with the investigation of convergence, showing that the LDE together with 

the resummation technique can exactly reproduce the large- iV result already at the first non-trivial order, provided 
that it is applied to a specific approximation fully exploiting the infrared limit properties. Having explicitly shown 
the LDE convergence properties within a limit where an exact result exists we turn our attention, in Sec. V, to the 
realistic finite N case where no similar infrared a ppr oximation is available. Using the same resummation method, the 
order-(5^ results for ci and C2 obtained in Ref. |l8l | are improved and the results obtained seem to converge to the 
lattice Monte Carlo estimates of Ref. 14j. Our conclusions are presented in Sec. VI. For completeness and comparison 
purposes we also include two appendixes, one where the original large- iV derivation is reviewed and another one to 
detail useful properties of the large order behavior of the LDE. 



Let us start our work by reviewing the application of the LDE method to our problem. The LDE was conceived 
to treat nonperturbative physics while staying within the familiar calculation framework provided by perturbation 
theory. In practice, this can be achieved as follows. Starting from an action S one performs the following interpolation 



where Sq is the soluble quadratic action, added by an (optimizable) mass term 77, and S is an arbitrary parameter. 
The above modification of the original action somewhat reminds the usual trick consisting of adding and subtracting 
a mass term to the original action. One can readily see that at 5 = 1 the original theory is retrieved, so that S actually 
works just as a bookkeeping parameter. The important modification is encoded in the field dependent quadratic term 
Sq^t]) that, for dimensional reasons, must include terms with mass dimensions (77). In principle, one is free to choose 
these mass terms and within the Hartree approximation they are replaced by a direct (or tadpole) type of self-energy 
before one performs any calculation. In the LDE they are taken as being completely arbitrary mass parameters, 
which will be fixed at the very end of a particular evaluation by an optimization method. One then formally pretends 
that S labels interactions so that Sq is absorbed in the propagator whereas 6S0 is regarded as a quadratic interaction. 
So, one sees that the physical essence of the method is the traditional dressing of the propagator to be used in the 
evaluation of physical quantities very much like in the Hartree case. What is different between the two methods is 
that, within the LDE the propagator is completely arbitrary, constrained only to cope with the so-called direct terms 
{i.e. tadpoles) within the Hartree approximation. So, within the latter approximation the relevant contributions are 
selected according to their topology from the start. 

Within the LDE one calculates in powers of S as if it was a small parameter. In this respect the LDE resembles 
the large- iV calculation since both methods use a bookkeeping parameter which is not a physical parameter like the 
original coupling constants and within each method one performs the calculations formally working as if — > 00 
or (5 — > 0, respectively. Finally, in both cases the bookkeeping parameters are set to their original values at the end 
which, in our case, means S — I. However, quantities evaluated at any finite LDE order from the dressed propagator 
will depend explicitly on 77, unless one could perform a calculation to all orders. Up to this stage the results remain 
strictly perturbative and very similar to the ones which would be obtained via a true perturbative calculation. It is 
now that the freedom in fixing r] generates nonperturbative results. Since 77 does not belong to the original theory 
one may require that a physical quantity ^^'^^ calculated perturbatively to order-J*^ be evaluated at the point where 
it is less sensitive to this parameter. This criterion, known as the Principle of Minimal Sensitivity (PMS), translates 
into the variational relation Isjl 



II. LDE AND THE INTERPOLATED EFFECTIVE SCALAR THEORY FOR BEC 



Ss = SS+{l-S)So{v) , 



(2.1) 



= . 



(2.2) 
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The optimum value fj which satisfies Eq. 12. 2() must be a function of the original parameters including the cou- 
plings, which generates the nonperturbative results. Another optimization procedure, known as Fastest Apparent 
Convergence criterion (FAC) (see also Ref. [S^), may also be employed. It requires, from the fc-th coefficient of the 
perturbative expansion 



$(fc) ^ 



i=0 



(2.3) 



that 



(5=1 







(2.4) 



1) in Eq. (|2.3() equal to zero. For the interested reader 
l33l | provide an extensive (but far from complete) list of 



which is just equivalent to taking the fc-th coefficient (at 6 = 
Refs. [H III El IIS El 111 111 111 111 117, 28, 2l|Ml^ 

successful applications of the method to different problems. 

It is important to recall that the basic reason for the convergence of the LDE method in the quantum mechanics 
case (anharmonic oscillator energy levels typically) '25', '2I HI relies on the fact that the LDE modifies perturba- 
tive expansions in such a way that the PMS or FAC optimized values of the initially arbitrary mass parameter (the 
equivalent of ry here) essentially follow, at large perturbative orders, a pattern of rescaling this mass with the pertur- 
bative order, which is such as to compensate the generic factorial growth |35| of the original perturbative expansion 
coefficients at large orders. As we will see in Sec III, in the present BEC case the relevant perturbative expansions 
do not exhibit factorially growing coefficients, but nevertheless the reasons for convergence of the LDE share some 
similarities with the latter cases, since the LDE followed by PMS also introduces at sufficiently large order a certain 
scaling behavior with the perturbative order, in such a way as to modify (extend) the convergence radius of the 
original perturbative series. 

Let us now write the interpolated version of the effective model described by Eq. (|1.8|) . Before doing that let 



us rewrite rba 



A where ^ is a mass counterterm coefficient needed to remove the UV divergence from 



the self-energy. This counterterm is the only one effectively needed within the modified minimal subtraction (MS) 
renormalization scheme adopted here for the evaluation of the relevant Feynman diagrams contributing to the self- 
energy, as explicitly shown in Ref. Then, one can choose 



5o = 2 [m' 



obtaining 



Ss 



|V0r + -ry"r + -(r-r;^)0^ + 



5u 



(2.5) 



(2.6) 



where As represents the renormalization mass counterterm for the interpolated theory, which depends on the param- 
eters 1] and 5. It is important to note that by introducing only extra mass terms in the original theory the LDE does 
not alter the polynomial structure and, hence, the renormalizability of a quantum field theory. In practice, the original 
counterterms change in an almost trivial way so as to absorb the new 77 and 5 dependence. The compatibility of the 
LDE with the renormalization program has been shown in the framework of the 0{N) scalar field theory at finite 
temperatures, in the first work of Ref. plj . showing that it consistently takes into account anomalous dimensions in 
the critical regime. Note also that we have treated r as an interaction, since this quantity has a critical value (rc) 
that is at least of order 5. 

Requiring the original system to exhibit infinite correlation length at the critical temperature, means that, at T,. 
and 5 = \ (the original theory), the full propagator G^^"'{p), given by 



must satisfy G^^\Qi)~^ = 0. This requirement implies 



(2.7) 



<5r(^') = -SW (0) , 

which is equivalent to the Hugenholtz-Pines theorem applied to the LDE. 



(2.8) 
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III. LDE EVALUATION OF ((/>2)i*^ IN THE LARGE-iV LIMIT 



Let us now turn our attention to the explicit LDE evaluation of (0^)ti in the large-A^ limit. In practice, the 
large- A'^ evaluation can be performed in different fashions which include the conventional order by order perturbative 
evaluation and the more economical closed form evaluation in which the whole large-iV series is resummed. The 
first, purely perturbative method in the standard Feynman graph way is also the only possible one concerning the 
finite N calculations, where different classes of diagrams contribute. The second technique is usually employed in 
approximations such as Hartree and 1 /N where it is possible to sum a certain class of graphs based on the type of loop 
terms they contain. Having resummed a given class, one may easily obtain a perturbative result by expanding the 
series to a given order in Su. Of course, both methods must lead to equivalent analytical results but, as we shall see, 
the final numerical results can be different at high perturbative orders. This is due to the fact that both perturbative 
expansions contain coefficients numerically produced. Since our optimization procedures may be sensitive to the 
numerical precision of those coefficients it will be instructive to compare them in detail. Finally, one can move one 
step further by obtaining a series with exact coefficients that allows for a fully analytical investigation. This is made 
possible by considering an approximation which avoids complicated integrals appearing in the exact calculation due 
to the presence of dressed propagators in terms of self-energies, which are usually cumbersome beyond some given 
order. Such a simpler series, with exact coefficients, is obtained if one considers typically the physically motivated 
deep infrared behavior of the dressed scalar propagators. In this section we explore these three possible evaluations. 



A. Standard Perturbative Evaluation of {(j>^)^u^ 

Let us evaluate {4>'^)u^ in the usual perturbative way. The relevant contributions, in the large- limit, are shown 
in Fig. 1. Using the full propagator one may write this quantity, at the critical point, as 



N 



(2^)3p2 + (^*)2 



p2 ^ (^*)2 



(3.1) 



where r/* = r/\/T~~J. Note that with this prescription one only has to evaluate diagrams which would appear in a 
usual perturbative calculation since the quadratic Si]^ vertex is automatically taken into account when rj* is expanded 
to the relevant order in 6. 

One can express the large- iV calculation more conveniently, in the generalization of Eq. I|1.8|) to 0{N) symmetry, 
by considering u — u' /N. In this case the nontrivial contributions, in the large- iV limit and expanded to LDE order 
fc, are given by 



-N 



(Pp [E(")(p)-I](")(0)] 



(27r)3 p2 + (^*)2 J (27r)3 [^2 + (^*)2]2 

where we have used Eq. (|2.8|) and E^"^ denotes the n-bubble self-energy given by 



(3.2) 



N 



5u' 



6 



n+l 



dH 



(2^)3 (/2 + {r,* f ) 



dh 



(2^)3 (s2 + (r;*)2) [is + p - ly + i-n*)^^ 



(3.3) 



which is then of order k — n + 1 in S. Note that the mass counterterm is a redundant quantity in the evaluation of 
{4>'^)'u^ because this quantity depends on the difference 



i:[:i{p) - i](^^)(o) = [i](^(p) + s(^^(p)] - [s|;'^(o) + Y.'ii^m , (3.4) 

where I]|["^(p) is the divergent self-energy. For a general renormalizable theory, the quantity Yi^^l\p) represents all 
counterterms associated with the parameters of the theory (such as masses and coupling constants) as well as the wave- 
function counterterm associated with any eventual momentum dependent pole. At the same time, E^"'' (0) involves 
the same counterterms except for the wave-function one. However, as we have already emphasized previously, in the 
three-dimensional case the only type of primitive divergence requires only a mass counterterm, which is the same for 
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^d?v(p) ^cUv ('-')• This means that in our case, S^"^(p) — E^"^(0) is always a finite quantity. It turns out that this 
quantity is also scale independent as discussed in Ref. '18^ . Since however the individual contribution S^"^ (p) contain 
a divergence, we regularize all diagrams with dimensional regularization in arbitrary dimensions c? = 3 — 2e, where in 
the modified minimal subtraction (MS) renormalization scheme, the momentum integrals can be written as 



(27r)3 



47r 



(27r)" 



(3.5) 



where M is an arbitrary mass scale and je — 0.5772 is the Euler-Mascheroni constant. 

Then, from the use of standard Feynman parameters for the integrals over momenta, we can write the general form 
for each of the two terms in Eq. (|3.2I) . that depend on the self-energy. The first of such term can be expressed 

as 



N 



d^p I](")(p) _ (-5u')"+^r[n(l/2 + e) + 2e] /e''M2y^"+^^ 



(27r)3 [p2 + (ry*)2]2 3(6?7*)" (47r)3/2(n+2) 

Jo Jq 
X /'dai[ai(l-ai)]-(i/2+e)... T da„[a„(l - a„)]-(i/2+e) 

jQ Jo 



Jo 



(3.6) 



where 



^(x,7,«.,/5,) = |x + :^^^^3^ 



7+(l_,)( Ji^+Mi^ 

ai(l - ai) 02(1 - Q!2) 



{f3iP2 ■ ■ ■ Pn^l) 

an(l - an) 



-[n(l/2+c)+2£] 



At the same time, the p = term is given by 



(3.7) 



TV 



d^p I](")(0) (-5w')"+^ r[n(l/2 + e) +e- 1/2] /e^M2y("+^^ 



(27r)3 [p2 + (r,*)2 



3(6??*)" 87r(47r)3/2(»+i) \^ (^*) 



X / d7(l-7)"(i/2+')-i[7(l-7)]i/2— "(1/2+.) 

dai[ai(l - ai)]-(i/2+.) da„[a„(l - a„)]-(i/2+') 

Jo 

X / dp^pr'mi - f d/32/5r'[/3i/32(i - m-'^^ . 



X / d/?„_i[/3i/?2.../3„-i(l-A.-i)]"'/'e(7,a»,/3j) 
Jo 



(3.8) 



where 



7(1 - 7) 



7 +(1-7) 



(l-/3i) , Piil-M 



ai{l — ai) 0:2(1 — 02) 



(/3l/32.../3n-l 

a„(l - an) 



-[n(l/2+e)+e-l/2] 



(3.9) 

It is not very difficult to see, by counting the superficial degree of divergence in Eq. H3.3|l . that the only ultraviolet 
divergence shows up in the one-bubble {n = 1) contribution. In Eq. H3.6|l the UV divergence for this case hides in the 
term x(l — T^)-3/2-i-e+ra(i/2+e) g^j-^j appears explicitly upon integration by parts over x- After that one can take the 
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usual expansion in powers of e and perform a numerical integration over the Feynman parameters to obtain for the 
first, non-vanishing, term in Eq. (|3.(jl) the result 



p 



i + 61nf — I +2-41n2 



[p^ + {-q*fY V* 18(87r 

In the p — case the pole shows up in the gamma function which becomes r(2e) for n = 1. Integration yields 



(3.10) 



[p2 + (,y*)2]2 ^, 18(87r)3 



7+61n(^)+2 + 41n(2/3) 



(3.11) 



The last two equations also reproduce the results found analytically in Refs. [T7','36'|. As already mentioned, although 
Eq. (|3.1l)|l and Eq. (|3.11|l diverge, their sum is finite and scale independent. Together, they give the contribution 

Jp + V* 18(8^)3 

All higher loop contributions are finite and one can safely take e = in Eqs. (|3.6(l and l|3.8(l . Note that in the above 
perturbative series, that is generated order by order from the standard Feynman graph procedure, the first non-trivial 
contributions start at order S'^u'^, due to the fact that the first order expansion term, linear in Su {i.e. n = in 
Eq. 1)3. 3|l ). is automatically canceled as a consequence of the Hugenholtz-Pines theorem, Eq. (|2.8(l . The number of 
Feynman variables at each order is fc -I- 2 for Eq. H3.6|l and fc + 1 for Eq. (|3.8|) . We then get the order-J^" result in 
the large- iV limit 



a2\(20) _ 



Except for the first coefficient, where we have the exact result from Eq. H3.12|l . all the other coefficients for i > 2 
can only be obtained numerically. One well known numerical routine that can be used to evaluate the i-dimensional 
integrals over the Feynman parameters in Eqs. (|3.t)|) and H3.8f) is the Monte Carlo multidimensional integration 
routine VEGAS [s^]. However, one should bear in mind that VEGAS may not be so reliable for a very large number 
of dimensions, since VEGAS, as a Monte-Carlo integration method, inherently makes use of finite numbers of points 
and iterations and these cannot be increased indefinitely in practice in order to improve precision. So for integrals 
with a very large number of dimensions (for example the last coefficient Cig involve a 39-dimensional integral) the 
VEGAS may lead to wrong estimates for both the numerical value of the integral and the corresponding error bar 
estimate obtained from the code (which also depends on the number of points and iterations used). Fortunately, as 
we will see below, all the terms contributing at large N can be computed alternatively in a much easier way (and to 
arbitrary precision) in terms of one-dimensional integrals, thus assuring a much better precision for the results and 
we do not need to worry about any specific detail of any numerical routine to evaluate Eqs. H3.6|l and (|3.8|) . The 
coefficients obtained this way, which we denote by Ji and are given in the following Subsection, will be the results 
used in all of our large- iV calculations. These Ji coefficients can be contrasted then to the results obtained from 
e.g. with VEGAS, that we show here for illustrative purposes only, obtained using 10^ points with 100 iterations 
fixed VEGAS parameters: Ci = (7.249 ± 0.001) x 10-^ C2 = (2.050 ± 0.003) x lO^^, C3 = (6.32 ± 0.01) x 10-^ 
C4 = (2.048 ±0.003) X 10"^ C5 = (6.85 ±0.01) x 10"", Cg = (1.709 ± 0.002) x IQ-^^, C7 = (3.561 ± 0.006) x 10"", 
Cs = (6.48 ±0.01) X 10-l^ Cg = (1.063 ±0.002) x 10"", do = (1.560 ±0.005) x lO'^^ Cn = (2.59 ±0.04) x lO'^i, 
C12 = (5.00±0.09) X 10-23^ Ci3 = (9±2) x 10-^6, = (5±l)x lO'^s^ Cig = (5±1) x 10-3°, Cie = (1±7) x 10-33, 
Ci7 = (2.3 ± 0.6) X 10-34, C18 = (2 ± 2) X 10-3^ and C19 = (5 ± 5) x 10-39. Note that starting with C13 the errors 
increase considerably as the dimension increases, as expected with a fixed number of VEGAS parameters. 



B. Closed Form Evaluation of 

Let us now write the whole large- iV perturbative series in a closed form which resembles the usual 1 /TV resummation 
of Feynman graphs. For completeness and comparison purposes we re-derive, in Appendix A, the original large-iV 
result found by Baym, Blaizot and Zinn- Justin (IQj . Note, as already emphasized, that one basic difference between 
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the original large- A'' calculation and the LDE one is that the latter automatically introduces an infrared regulated 
propagator, from the explicit mass term 77. Apart from its main purpose of defining in that way the relevant LDE 
series in 5u/rj, cf. Eq. H3.13|l . this also has the advantage of explicitl y re gularizing the intrinsic infrared divergence 
of the corresponding expression of Tc in the original calculations "lOl lll|. However, similarly with the latter, there 
still remain some subtleties with this closed (resummed) form of the perturbation series, related to the fact that the 
integrals over momenta are not absolutely (UV) convergent, as we shall examine below. Thus, after applying the 
Hugenholtz-Pines theorem, and summing all the leading large- contributions shown in Fig. 1, one obtains for the 
expression equivalent to Eq. 1)3. 2|l the result 



SuN f (fp (fk 



3 J (27r)3 (27r)3 [p2 + (,y*)212 



(fc +p)2 + (ry*)2 /j2 _^ (^*)2 



where 



jfq 1 

(2^ [g2 + (77*)2] [(k + qf + {ri* 



1 / k 

— - arctan — 
47rfc \2r]* 



, (3.14) 



(3.15) 



with k = |k| and similarly for p, q in Euclidean space. Contrary to the corresponding expression in the massless case 
rj = (see Appendix A) , here it appears not possible to integrate exactly Eq. H3.14|l due to the non-trivial dependence 
in k and 77* of the resummed propagator B{k,ri*). But one can at least still do exactly the first integral over the 
momentum p. In Eq. (|3.14fl . the integral over p is finite in ci = 3, and can be easily performed to give 



/,2\(«) ^ _ 



Nt]* SuN 1 



An 



(fk 



3 (87r?7*) J (27r) 



6uN 
~6~ 



n -1 r 



B{k,r^*) 



fc2+4(77*)2 fc2 + (r;*)2 



(3.16) 



while the remaining fc-integral can be performed numerically. After a little algebra one gets 



47r 

where the Ji coefficients are given by 



i—1 ^ ^ 



167r3 



dz 



(z2 + l)(z2+4) 



(3.18) 



with 



A{z) — - arctan — , (3.19) 

and z = k/rf. Analytically, the two ways we have presented for obtaining the perturbative evaluation of {<jP')u^ are 
expected to be equivalent and any difference may only arise from the numerical evaluation of the Ci and Ji coefficients. 
In this respect one expects Ji, which are evaluated from one dimensional integrals, to be more accurate than d and 
it will be instructive to compare both results. We have numerically evaluated Ji with both MATHEMATICA 's^l 
and Maple, where we can compute the integrals with arbitrary precision, using diverse integration routines available 
in both, in order to check the reliability of the results. The first nineteen values obtained are Ji = 7.24858 x 10"^, 
J2 = 2.04919 X 10-^ J3 = 6.32139 x 10-^ J4 = 2.04829 x 10-^ J5 = 6.85295 x 10"", Jg = 2.3454 x 10-^2^ 
J7 = 8.16524 X IQ-^^, Js = 2.88069 x IQ-^^, Jg 1.02726 x IQ-^^, Jio = 3.69589 x IQ-^^, Jn = 1.33955 x IQ-^^, 
J12 = 4.88611 X 10-21, = 1.79203 x IO-22, j^^ ^ 6.60406 x 10-24, J^g = 2.44405 x 10^25, j^g = 9.07903 x 10^27, 
Ji7 = 3.38400 X 10-28, = 1.26514 x lO-^^ and J19 = 4.74288 x 10-^^ Comparing d to J, one sees that the 
multi-dimensional VEGAS routine produces accurate results up to i = 5 but the values quickly deteriorate at large 
orders. This is a consequence that the VEGAS routine does not handle well integrals with a too large number of 
dimensions, for relatively (and computationally viable) small number of points and iterations, as explained at the 



11 



end of Sec. III. A. Therefore, in this work, we shall consider only the series with the more accurate Ji coefficients, as 
computed from Eq. I|3.18|l . 

However, one notes that Eq. H3.17(l displays a more significant difference with respect to the "direct" perturbative 
series in Eq. 13.13|l . which is due to the presence of the extra first order term in (5u, independent of 77. Remark that 
this contribution is just the opposite, in sign, of the exact large- iV result ll^] (see Eq. IjAllI) in Appendix A). This 
apparent difference between the two approaches to the large- A'^ perturbative series deserves a detailed discussion to 
which we now turn our attention on. In fact, the above difference is only a consequence of integrating expression 
(|3.14|) over p first. Namely, if performing the expansion in powers of 5, and then integrating first over k in Eq. (|3.14|l 
(which is formally equivalent to what is done in the standard perturbative order-by-order graphical approach in the 
previous subsection), the linear term in 5u automatically cancels out. More precisely one obtains in this case an 
integral of the type (^ just denotes an arbitrary mass parameter here, that is equal to r]* in the LDE calculation) 



and in dimensional regularization the two terms in the bracket just cancel out, as can be seen by making a shift 
k ^ k — p in the first term. (One can also check with a standard cutoff regularization that the integral in Eq. (|3.2(J|) 
gives a zero result, though it is a less immediate calculation). Therefore depending on the order in which the two 
integrals in Eq. 13.14|l are performed one may get different results, which is precisely the manifestation of an ambiguity 
due to the fact that the integrals are not absolutely (UV) convergent, as pointed out in Refs. Actually, this 

problem is more basically rooted in the fact that in obtaining Eq. (|3.14l) one has formally resummed a series containing 
UV divergences, considering, e.g., the separate contributions in the last parenthesis of Eqs. H3.14|l (see also Eqs. (|Aip 
and HA2|I in Appendix A). Therefore, the actual point is that one is not allowed (in principle) to exchange the 
perturbative, all order summation, with integration, which in our calculation is reflected in the different resulting 
perturbative series. Going back to the standard perturbative expansion, as performed in the previous subsection 
(Eqs. (|3.2(l - (|3.8|l ) the perturbative parameter expansion in powers of 5 is made first, and the UV divergence (which 
only appears at first non-trivial order as discussed there) may be taken care of by the standard renormalization. 
On the other hand, if we formally re-expand Eq. H3.14|l in a power series in u, we can immediately see that the integral 
defining the coefficient of the ffist order term linear in u5, originating from [1 + (5uN /Q)B{k, ?7*)]~^ ^ 1, has momenta 
routing that cannot be consistent with the actual perturbative graph: rather, considering for instance the first term 
of the last parenthesis in Eq. H3.14|l , one should have p + fc — > p for consistency (see Fig. , since this first order term 
in 5u implies that the resummed propagator (dotted line) is pinched to a point, so that there is no fc-momentum flow. 
This would give again a zero result for the coefficient of u, just for the same reasons as Eq. (|3.2U|) is vanishing, while 
formally performing the integral with p + k instead gives the A^u/ (967r^ ) term 2. To summarize, while it will appear 
very convenient to formally perform the integral over p first in Eq. (|3.14(l , it is only consistent provided one subtracts 
the spurious linear term from the naive result, Eq. H3.17|) . The correct perturbative series thus reads 



+ ' (3.21) 

i—l ^ ^ 

which has the same form as Eq. (|3.13|) . An expression similar to Eq. H3.21|l was also found by Braaten and Radescu 
inRef. I^. 



C. Asymptotic infrared and ultraviolet behavior of {cf>^)u ^ 

Before considering the relevant BEG perturbation series, Eq. (|3.13|) . or equivalently Eq. H3.21|) . in the large- A^ limit, 
let us recall some expected general properties of the large order perturbative expansions, as seen from a diagrammatic 
point of view. This digression will emphasize an important difference between the generally expected large order 
behavior of perturbative series in most renormalizable models and the behavior of the above BEG specific series. 



^ Remark that in the original calculation lij corresponding to take = in Eq. 13.141 . this ambiguity problem has been consistently 
solved simply by using dimensional regularization (see Appendix A for details): in contrast, the rather subtle point is that when rj 
in Eq. I3.14i it appears at first perfectly consistent to start with the p integration. Nevertheless, this does not give the correct result, 
independently of whether one uses dimensional regularization or not, until one correctly identifies what the actual perturbative series 
in 5u/r] should be, as explained above. 
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In field theory one has to face the problem of the perturbative series being often only an asymptotic (non-convergent) 
series which, in most cases, is due to factorial growing perturbative coefficients at large orders. If the coefficients 
are of the same sign, order by order, those series are not even Borel summable j89t I4(tI |. In practice, this means 
that the perturbative expansion alone does not define uniquely the physical quantities being expanded, so that the 
series has to be complemented by intrinsically nonperturbative contributions, containing typically terms with an 
exponential dependence in the (inverse) expansion parameter M- This is problematic because, apart from the 
special case of exactly solvable/integrable models, in most theories those nonperturbative terms are at best known 
only on phenomenological grounds. However, to investigate the large order behavior of the perturbative series (and 
therefore guess at least the form of nonperturbative missing contributions), it is often sufficient to consider a class 
of approximated graphs, expected (and proved in some specific models) to give the dominant contributions to the 
perturbative coefficients at large order. In dimension c? = 2, 3 and 4 renormalizable theories, such dominant graphs are 
typically given by the next-to- leading term in a 1/A'^ expansion, where, roughly speaking, the matter fields are in a N- 
vector representation, and it is sufficient to consider the asymptotic behavior in the bubble (1/-/V) approximation of the 
relevant Green functions. For instance, in a theory with a renormalized coupling one obtains, after renormalization, 
a scale-dependent running coupling. A very sketchy Green's function calculation in the above approximation involves 
(after renormalization) typical momentum integrals of the form 



-1 

(3.22) 

where F{q^) is model-dependent and characteristic of the Green's function considered, g{M) is the running coupling, 
/So is the first order renormalization group (RG) /3-function coefficient: dg{M)/dhiM = Pog'^{M) + ■ ■ and M is an 
arbitrary renormalization scale. When formally expanding Eq. H3.22|l in a perturbative series in g, one gets integrals 
of the form 



dqF{q^)g{M) 



/3o.g(M)ln^ 



E 



>+i)(M)(/?o)^ r dq^F{q^)\n^(^] , (3.23) 



which leads to a factorial behavior p\ at large p. More precisely, F{q'^) can be expanded in a power series in q^ (l/?^) 
in the infrared (respectively in the ultraviolet), so that Eq. (|3.23|l gives series of the form ^3 ~ g^^'^^H~Poy p! 
or ^ 3^^^^-* (/^o)^ for large p, respectively, for the infrared and ultraviolet behavior. Considering for example an 
asymptotically free theory, i.e. with /Jq < 0, one obtains a non sign-alternated series thus non-Borel summable, as far 
as the infrared behavior is concerned. This fact reflects the important infrared sensitivity of such theories. 

Now, a drastic difference between the previous illustration of theories with a renormalized coupling and the effective 
BEG 0^ model in three dimensions considered here, is that for the latter only the mass is renormalized, so that the 
coupling is finite and dimensionful, as pointed out previously. From this, and following the above line of reasoning, 
one expects that the relevant BEG perturbative series coefficients at large orders do not display any factorial 
behavior. Therefore, a more convergent series should appear, as is confirmed for instance by the form of the exact 
large- iV perturbative series in Eqs. H3.13|l . (|3.21|) . whose coefficients appear clearly not very different from those of 
a geometric series. Also, from the above general considerations, a similar behavior of the series is expected as well 
for arbitrary N. Thus an interesting question is whether one could obtain from such large order behavior estimates 
a sensible approximation of the exact series in Su/rj, that would be relevant within the LDE method. 

Let us therefore investigate some analytically simpler but physically motivated approximations (expected to be 
asymptotically dominant) of the large order behavior of the power series in Su/rj, as generated from Eq. H3.16|l . 
B{k,r]) behaves as 



B(fc,77*)~<j (3.24) 

for the IR (fc <C ??*) and UV (fc 3> r/*) limits, respectively. This means that the deep IR behavior of the series should 
be essentially given by a fc-independent term, [1 + 6uN / (ASmf)]'^ , replacing the corresponding bracket in Eq. H3.14|l . 
Retaining only this simplest IR behavior, the remaining integral over k becomes straightforward and accordingly we 
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obtain^ the relevant IR approximated series as a simple geometric series 



Nrj* SuN 



Att 247r77* 



SuN 



2t]* 
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SuN 
967r2 



Ntj* SuN 



SuN 

AS-KT]* 



(3.25) 



where the first order term, independent of rj, has been explicitly subtracted from the naive integral result Eq. (|3.16() 
for consistency, as discussed in detail in the previous Subsec. III.B (note that the purpose of the last parenthesis in 
the first equality in Eq. (|3.25|) is to retain, for clarity, the separate contributions of the two propagator terms in Eq. 
(|3.16l) ). Note also that Eq. (|3.16l) is UV finite, so that the result Eq. H3.25|l is independent of the integration method 
used, so that either dimensional regularization or another integration method leads to the same result. 

Similarly, we can still integrate Eq. H3.14|l exactly when taking the UV limit of the propagator, B{k,jf) (8fc)^^ 
from Eq. H3.24|l . One obtains 



[-J/(7 + 4y^)-81n(2)(l + y^)-61ny] (l + 5,^ + 4^^)-^ " (^ ' ^^'^^^ 

where y = ASrj* /{N Su). When performing the p integration first, all integrals are UV finite in d = 3 and again the 
r] ^ mass explicitly regularizes the IR divergences. Note for instance that in both, Eqs. (|3.25() and H3.26|l . the l/rj* 
and J]* from the first and second integrals respectively, cancel out: this is completely analogous to the cancellation in 
the original calculation lOj, in dimensional regularization with no IR cutoff mass, of the pole l/((i — 3) (see Appendix 
A for details). 

Next, remark that by taking the r] limit of expression H3.25|l or H3.26|l . one recovers the correct 1/A^ exact 
result Eq. HAllI) . In other words, we see here that the massless limit 77 of the LDE series consistently reproduces 
the exact large- iV result, which is expected since the LDE for S 1 plays the role of an infrared regulator. But 
this check is important as regards the question of the possible convergence of the LDE series to the exact result, 
once a non-trivial optimization of the LDE series, with respect to the arbitrary remaining mass parameter 77, will be 
performed, to be examined in the next section. 

With this aim, it is instructive to re-expand, in a power series of Su/rj*, the above two different IR and UV 
approximations for A(<jP)'^^\ First, taking the IR limit Eq. (|3.25|) gives the geometric series in Su/rj*: 



(*^>s^-4f-^|:«(-|^)" (3-) 

where Gi = [(647r2)(87r)*]-i. Numerically, the first five G^ coefficients are d = 6.2991 x 10~^ G2 = 2.5063 x 10"^ 
G3 = 9.9724 X 10"®, Gi = 3.9679 x 1 0"^ a nd G5 = 1.5788 x 10-^°, which are interesting to compare with the 
corresponding exact coefficients in Eq. H3.21(l . One can see that the first low order coefficients of Eqs. (|3.27(l and 
(|3.21l) are of very similar magnitudes, and we have further checked that significant departures {i.e. about an order of 
magnitude or more) between the IR-approximated and exact large- iV perturbative coefficients only occur at rather 
large, greater than 0{S^^), orders. In other words, one expects from this comparison the IR-approximated series, 
which has a convenient and simpler geometric form, to be a very good approximation of the exact large-iV series. 
This is a strong indication that the detailed non-asymptotic (infrared) behavior of the scalar propagator should play 
essentially no role, as could be physically expected on general grounds, and as will be fully confirmed by our numerical 
investigation below. 

The same expansion for the UV limit H3.26(l reads similarly 



UV 



Nrj* S uN 
'^tT ~^ 3 



-1.7465 X 10~ 



Su 
3^ 



2.4737 X 10"^ 



(3.28) 



where one can see in contrast that the coefficients are already quite different from the exact series Eq. (|3.21l) at low 
order, so that we may expect the asymptotic UV limit of the propagator to give a less sensible approximation than 



^ Note that, in close analogy with Eqs. i3.22i . i3.23i . we only approximate the non-trivial resummed propagator B{k,ri*) according to 
Eq. 13.241 . and keep the exact k,r]* dependence of the remaining integrand. 
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the IR one. For completeness it is useful to consider alternatively the direct evaluation of Eq. H3.14|l {i.e. not as a 
perturbation series in Su/r]). Taking thus the exact expression for B{k, rj*) instead of its simpler IR or UV limits, the 
^-integration can be performed only numerically (note that it is still IR and UV finite in d = 3). Those (numerical) 
results for the reference value u — 1 and as function of 77* are illustrated in Fig. [3 As one can sec, both the IR 
and UV approximations, as given in Eq. (|3.24() . have a behavior that is somewhat different from the exact function 
Eq. H3.14|l for very small 77*, although all expressions correctly give the exact result at rj* — 0. However, though it is 
not visible on Fig. 2, the IR and UV approximations appear to be very good approximations of the exact function for 
larger 77* 0{1), that is in the range where their respective perturbative expansion forms start to be valid. 

To summarize this subsection, introducing an IR regulator mass 7^ as is done from the LDE procedure, together 
with the deep IR limit of the propagator, Eq. 13.24|l in Eq. I|3.14|l leads to perturbative series that are very close to the 
exact large- TV one in (|3.21|l . The only subtlety when implementing the LDE within the convenient resummed large- A'^ 
closed form Eq. (|3.14|) , is to remind that one should be careful in exchanging the perturbative series summation with 
integration, since the resulting integral, Eq. (|3.14() . is not (absolutely) UV convergent. The consequence is that, e.g., 
the 77* — > limit of expression (|3.16|l does not commute with taking the limit rj* ~> before performing any integration, 
i.e. as is done in the original large- iV calculation (see Appendix A). This may be considered as a reminiscence of the 
infrared sensitivity of the theory, even if it is not as severe as the above mentioned models with a running coupling, 
leading to divergent series with factorial growing perturbative coefficients. 

IV. LARGE-iV OPTIMIZATION, RESUMMATION AND CONVERGENCE PROPERTIES 

As discussed in the introduction, the study of LDE convergence properties in the BEC problem is much more 
complicated than in the pure anharmonic oscillator case 2^, 113 Ell • In principle, both models are described by 
a scalar model in one and three dimensions but in the BEC case, the model is used to study a phase transition. 
However, if the LDE works, one expects that reasonable numerical results should be obtained, converging to the 
"exact" large-A result, ci — 87r/[3C(3/2)^/'^] ~ 2.328, evaluated in Ref. Before proceeding, it is useful to 

point out an essential aspect and potential difficulty of the purely perturbative expansions, Eq. (|3.21|) . or similarly 
Eq. H3.27|l . Note that these are inverse (alternated) series in the mass parameter rj, thus with an a priory finite 
convergence radius 77" -'^ (i.e. for u — 1 these series are (absolutely) convergent only for I77I > 77c) when considered as 
series for complex values of the arbitrary mass parameter 77. But, as discussed above, ultimately the exact result is 
expected to be recovered for 7; — s- 0. This situation is not much different from the anharmonic oscillator case, where 
typically the energy levels have perturbative expansions in powers of X/m^ (A being the coupling and m the mass) 
with moreover factorially growing coefhcients at large perturbative orders, but where nevertheless the LDE converges 
HE HE nil to the exact result, thanks to an appropriate rescaling of the mass parameter that is consistent with the 
PMS optimized solutions. We will examine here how the LDE procedure followed by the standard PMS optimization, 
Eq. H2.2|l . manage in fact to avoid this 77 — > potential problem with the basic perturbation series, which is one of the 
main results of the present paper. This is where the infrared approximation is a useful guide: while its perturbative 
form Eq. H3.27|) exhibits just the same feature as the exact series, Eq. (|3.21|) . the former geometric series is known 
to all orders, and obviously its 77 ^ limit is perfectly well-defined and gives the correct exact result, as can be seen 
by going back to its original form Eq. 13.25() discussed in Sec. III.C. But no such resummation is a priory known for 
the relevant non-trivial N = 2 series, where only the first few perturbative orders in u/77 are known and thus only the 
latter information can be used to define the LDE procedure. 

Thus, before considering the LDE of the actual BEC series, it appears very instructive to first examine the same 
LDE procedure performed on a simpler model which shares many similarities with the relevant BEC problem. We 
will see that the example below illustrates very well the basic reasons for the success (or eventually failure) of the 
LDE followed by the PMS optimization method in the general case, beyond the specific BEC problem considered in 
this paper. 



Let us examine the properties of the LDE and subsequent PMS optimization in a general case by considering the 
following function which admits a simple alternated geometric series expansion: 



A. A simple example of the LDE-PMS convergence 




(4.1) 
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which expansion form is almost Hke our IR geometric series Eq. (|3.27|) . for u = 1 and x = l/rj*, apart from overall 
different normalizations (compare e.g. with Eq. H3.21|l or (|3.27|) ). Clearly, the exact expression in Eq. 1)4. l|l tends to 
— 1 for X oo, thus the goal is to examine whether the LDE procedure followed by the standard PMS optimization, 
which only uses the series expansion form in Eq. H4.1|l at successive perturbative orders, is able to reach such a result 
and in which way. To make contact with the LDE series we consider the series in Eq. H4.1|l with x = uS/r]* (except in 
the —1/x term on the right-hand side of Eq. 1)4. l|l where we take x = l/?7*), followed by an expansion in power series 
of d in which one takes i5 ^ 1. The result of this LDE at arbitrary order k can be expressed entirely analytically in 
this case (see also Appendix B)"': 



/ 1 r- (-1)'' ^n T[l - n/2 + k] ,^ , 

'^(^ " "Z'^)' ' - -x^ ni/2 -kni + k] + 5/"'^ r[l- n + fc]r[l + n/2] ^'''^ 

and the PMS optimization performed order by order (up to order k — 10) is shown in Table I. From this table, we 
can draw several important observations: 



i) Applying the PMS optimization condition Eq. I|2.2|) to $(a;) generates complex fj solutions. In addition, all 

solutions can be arranged into families which span the complex plane. In general, a new family arises at even 
orders and this is signaled by a first member which lies in the real axis. This pattern was also found in the 
anharmonic oscillator applications '2^ and in the finite N applications to the BEC case 18]. 

ii) Despite the fact that the (naive) expansion at arbitrary finite order has a finite convergence radius < 1 (namely, 

the LDE at order k only uses the information of the right-hand side of Eq. (|4.1() , where the series is absolutely 
convergent for \x\ < 1), the LDE-PMS procedure clearly converges to the right result for a; ^ oo (equivalently 
ry — > 0) and the convergence is quite rapid in this case. 

iii) Though most of the PMS solutions do converge to the correct a; ^ oo result, there is a family of real solutions, 

SO, that clearly converges to a different result (~ 2). 

iv) Some of the (complex) solutions converge more rapidly than others. In particular, at order 3 one of the solutions 

is already very close to the exact result, and also one of the two real solutions that does converge to the exact 
result, which appears only at odd LDE orders, is not the best one. Note also that it is the real part which 
converges to the right result. However, to obtain those results one must consider all, real and complex, fj 
values. In fact, the imaginary parts of the complex optimized values fj play an essential role and for instance 
suppressing even a small imaginary part results in a completely wrong and unstable result for the optimized 
series. As we shall investigate in more detail below, this is because, for convergence of the series, it is the value 
of I a; I (respectively l??!^"'^) which is relevant. 

We shall now investigate in more details the basic reasons for the main LDE convergence results ii) above. First, 
note that in Eq. H4.1|l we added an extra term 1/a; ~ rj*, since it vanishes for a; — > cx). In the relevant BEC case, 
as already discussed in the previous Section, this term has a very clear physical interpretation, as it corresponds to 
the tadpole graph. In the present case it may be considered only as a mathematical trick, which plays an important 
role for convergence. Indeed, when suppressing this extra linear "mass" term, one finds for the LDE-PMS in place of 
the results shown in Table I: —0.5625 and — 0.6574 ± 0.153/ e.g. at LDE orders k = 2 and fc = 3 respectively; and 
at order fc = 10 the best solution is — 0.97208 ± 0.285/ while other solutions are still about 20% away from the exact 
result. Thus, the additional linear term clearly increases substantially the speed of convergence. First, it is clear that 
non-trivial PMS solutions already start at order 1, while they start only at order 2 when suppressing this term. But 
this is not the main reason for this faster convergence. Remark that at very large order, say k ~ 100, the numerical 
accuracy from both procedures {i.e. with or without the additional linear term) tend to become essentially equivalent 
(except that the incorrect solution ^ 2, see iii) above, is absent in the procedure with the linear term suppressed). 
However, it is evidently crucial to have a procedure such that the very first few LDE orders give already reasonably 



* Note a slight difference in Eq. 14.21 with respect to Eq. lB2t in Appendix B, namely the Gamma function ratio r[— n/2 + k]/r[—n + k] 
in place of r[l — n/2 + fc]/r[l — n -\- k] appearing in Eq. IB2I . due to the fact that Eq. 14.11 does not have the extra factor of u5 in 
front of the series. Because of this, non-trivial PMS solutions d<l?(x)('°' /dx = of Eq. 14.21 start already at LDE order fc = 1, see Table 
1. But apart from that, this difference only affects in very minor ways the qualitative behavior discussed below, in particular the large 
LDE order behavior of the simple series ^4.21 with respect to the corresponding actual BEC series. 
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reliable results, and also to gain an order in the PMS solution, when we will consider the relevant N = 2 BEC series 
for which only the first few perturbative orders are known. The essential role of the linear term is easily understood 
when expanding the (exact) result in Eq. H4.1|l for large x: 

such that for x ^ oo (77 ^ 0) the first order 0(1/0;) term cancels out. Now what happens is that this cancellation 
remains after introducing the LDE procedure, thus leading to a faster LDE-PMS convergence, see remarks ii) above, 
even though the LDE modifies the structure of the perturbative expansion and, as already mentioned, only uses the 
perturbative information from Eq. (|4.1|l . This can be understood on basic grounds when considering the large k 
behavior of expression H4.2fl : 



^(x = u/ri)'-'''^ ^ ——— + Y wo r , exp(a;^fc) Erfc(xk^/^) ~ 1 (4.4) 

fc->oo a;fc^/^r[l/2] r[l + n/2] j;fci/2r[l/2] v\ ) 1 \ ) \ J 

where Erfc{x) is the standard error function |4lj |. and this large order behavior in (|4.4|l is obtained after some 
algebra by using standard properties of the Gamma functions (see also Appendix B). Furthermore, it appears that 
this large k behavior of Eq. (|4.4|) exhibiting the scaling fc^/^ is rapidly reached: for instance the difference between 
(-l)'=/r[l/2 - k]/T[k + 1] in Eq. and k-^/^/n in Eq. (jOJ is already < 4% for fc > 3. Analyzing thus from 

Eq. 14. 4|) the large k behavior, one can derive by using the asymptotic expansion of Erfc{x) for the relevant limit 
X ^ 00: exp(a;^) Erfc{x) ~ 1/{^Jt:x) [1 + ©(l/x^)] (see also Eq. (|B5|I in Appendix B), that the first order term 
Oil/x/k^/"^) compensates exactly the linear term, giving in Eq. (|4.4|) <^{x) — 1 + 0(l/a;^/fc), while omitting the 
linear term gives instead: '^{x) '--^ — 1 + 0{\/x/k^/'^). Though the latter properties result from taking the somewhat 
extreme limit fc ^ 00, what is quite remarkable is that this behavior is already well observed at very low orders, as 
the above comparison of numerical PMS optimization in Table I illustrates. 

More generally, one can also understand from Eq. H4.4|l the main transformation operated by the LDE on the 
original series: while the latter had a finite convergence radius, for I77I > 1, the LDE provides an extra damping factor 
l/r[l + n/2] which allows to safely reach larger \x\ values (equivalently smaller \rj\ values) from the new perturbative 
series, so that eventually the relevant limit 77 — > may be approached. Indeed in practice, namely at finite k orders, 
all of the PMS solutions fj of Eq. H2.2|l corresponding to Table I (except the ones in the first column corresponding to 
the incorrect solution, as will be discussed below) tend to have smaller and smaller jryl values (though rather slowly 
decreasing) as fc is increased. 

A last remark on the LDE behavior of the simple example (|4.1|l concerns the occurrence of the incorrect PMS 
solution, as indicated in iii) above. Clearly, this results directly from the presence of the additional linear term: in 
absence of the latter, these extra solutions disappear from the LDE-PMS at any orders. Actually, these are reminiscent 
from the fact that Eq. (|4.1|l (before LDE is performed), has another extremum at a; = —1/2 [i.e. 77 = —2). In the 
most general case where one would have no a priory idea e.g. on the sign of the correct solution, this feature may be 
considered as a drawback of our procedure. But in fact, it is easy to get rid of this incorrect solution, simply because 
one knows that the solution we seek should be for ry — > 0. In contrast, the PMS optimized values at successive LDE 
orders fc corresponding to this extra solution are always such that \f]\is maximal, with respect to all other solutions. 
Furthermore, these |?7| values do not exhibit the expected trend towards smaller and smaller values: on the contrary, 
the corresponding \'q\ value is (rather slowly) increasing as fc is increased. 

The important point is that, as we will see below, all of the above properties will be exhibited similarly by the more 
complicated BEC series and are, therefore, a very useful guideline. In particular, the cancellation due to the additional 
linear term as above observed, the behavior with the LDE order fc of the PMS solutions, including the behavior of 
the incorrect PMS solution, all occur similarly in the more complicated cases of the actual BEC LDE series, with 
an expected much faster LDE convergence. These are important remarks concerning the LDE applications since, as 
already emphasized, the difference between our original applications Refs. '301 and the one performed in Ref. 
amounts to the fact that the latter authors optimize A(0^)('') in Eq. 1)1. 9|l which differs from (cfp)'-^^ by an equivalent 
linear term. 

In the next subsections, we switch to the study of the LDE performed on the actual BEC series in the large- 
limit. With this purpose we will consider ((/'■^)1'^'* as given by the closed form series, Eq. H3.21|l . and alternatively 
also its infrared approximation form, Eq. (|3.27() . These two quantities will be again optimized with the standard 
PMS criterion, Eq. (|2.2|) . For completeness we also show, with one example, how the alternative Fastest Apparent 
Convergence (FAC) method, Eq. (|2.4fl . generates similar results. We then present in Sec. IV. C a way to further 
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rearrange the LDE series, by fully exploiting the above explained behavior of the PMS optimized solutions, which 
eventually accelerates further the convergence due to the fact that it can recover directly the large LDE order behavior 
and consequent good 77 ^ properties of the infrared approximation, with the advantage of being applicable in a 
general case to more arbitrary series. This technique is then used to treat both the complete closed form and the 
infrared approximation of {(j)^)u \ as introduced respectively in Sees. III.B and III.C. For comparison purposes an 
alternative resummation procedure, based on Pade approximants |4^ . is also studied in Sec. IV. D. Finally, we analyze 
all results to conclude in Sec IV. E about the convergence structure of the LDE method. 



B. Standard optimization 

In order to perform a numerical analysis of the optimized solutions for (0^)1'^'' one first expands Eq. H3.13|l to the 
desired order in S. As usual, one then sets S = I before optimizing. This is the LDE part of the procedure. It can be 
performed either by explicit order by order expansion, or equivalently more formally following the general structure of 
the LDE expansion at arbitrary order k as presented in Appendix B, leading to the result Eq. ljB2p . Like the simple 

geometric series considered above, Eq. (|4.1f) . the equation for {<jP')'"u^ is basically a series in powers of {5uY / ■q^^~'^\ so 
that the PMS and Fastest Apparent Convergence procedures will also generate algebraic equations of order k whose 
mathematically acceptable roots form a set of optima f] values. 
Applying the PMS condition, Eq. to {(t>^)u\ 



drj 



= , (4.5) 



or the Fastest Apparent Convergence criterion, — (0^)1'^ = 0, computed dX 5 — 1 (which is again analogous 

of taking the fc-th coefficient in Eq. (|2.3(l . with <^^^^ = equal to zero) generates the optimal results. Table 

II lists, to order-(5^° all families obtained by applying the PMS optimization to the standard closed form expansion 
given by Eq. (|3.21|) . For completeness we compare, in Table III, results generated with the PMS optimization and 
the Fastest Apparent Convergence procedures for the same quantity. This table clearly shows that both optimization 
criteria generate compatible solutions that seem to converge to identical results. We advance that this pattern was 
also observed in all remaining cases considered in this work, so we shall consider only the PMS optimization, Eq. 
(|4.5|) . from now on. In addition, we can get more insight by considering our infrared approximated geometric series, 
Eq. H3.27|l . which has simple exact coefficients known to all orders. By applying the standard PMS, Eq. H4.5I) . to this 
simpler equation one obtains the results shown in Table IV. This table has exactly the same characteristics as Table 
II except for the positive real part numerical values, that are much closer to the exact large- A'' value ci ~ 2.32847. 

Let us try to examine at this stage the (eventual) convergence structure of the LDE as applied to the large- limit 
series, namely the results presented in Tables II-III and IV. In both cases the family structure is very similar to the 
one found for the simple geometric series example in the previous subsection, as well as in the studies that have proven 
the LDE convergence in the anharmonic oscillator For instance, we also see in the present case the appearance 
of a new family at even orders, such that all of them start with a real solution and become complex at the next order, 
the exception being FO which has only negative real solutions and FA which has only positive real solutions. As 
mentioned before, the complex parts arise as a consequence of solving the polynomial equations generated during the 
optimization procedure and are mathematically acceptable. We refrain from trying to attach any physical significance 
to the complex part of the optima ci values and instead of considering only completely real solutions we take the 
optimized Re(ci) as the relevant quantity for the evaluation of Tc, which is a strictly real quantity. Also, in the 
geometric series application above, we saw that only the real parts of the optimal solutions converge to the expected 
real value when x — *■ 00. Considering all complex solutions has also the advantage of giving a prediction at any order. 
Further, we note that all families whose real parts are positive start with values close to ~ 2.0 and then seem to follow 
similar patterns as the perturbative order increases. In contrast, considering FO would bring important qualitative 
changes regarding the critical temperature shift in relation to the ideal gas value, ATc = Tc — Tq, since the sign of this 
quantity has also been a source of controversy 9] for some time until recently, when most works started to predict 
positive values for the critical temperature shift. It appears clearly that FO is just the equivalent of the SO (wrong) 
solution discussed in the simplest geometric series case, which was eliminated because it does not correspond to the 
r] —f expected behavior. As we will see below, the same criteria allows to eliminate this solution without ambiguity. 
Pushing further to much higher orders we obtain for the exact 1/A^ series case Eq. (|3.21() . e.g. at LDE order 100, 
that most of the positive roots give a solution whose (real part) is ~ 2.2, as well as a few solutions which are far 
apart reasonable values, which we assume only reflect the numerical limitation of the problem at hand. Indeed, to 
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obtain those results the numerical procedure involves first to numerically solve the one dimensional integrals related 
to the Jj coefficients in Eq. 113.2111 . As emphasized, this has been done with great care with the maximum precision 
allowed by MATHEMATICA and/or Maple. Yet, one cannot expect the results to be completely stable once the 
LDE is carried out to very high orders. Though the integration in Eq. 13.18|l may be in principle done to arbitrary 
accuracy, the limitation comes about later in the process of finding the roots of high order polynomial equations, 
which is a notoriously unstable numerical problem in general. So it is rather the optimization procedure, as implied 
by the LDE/PMS which is to be carried out numerically, which appears sensitive to numerical accuracy at very high 
orders. Nevertheless from our analysis all results performed with MATHEMATICA appear under very good control, 
say until a LDE order of about ^ 50. It looks that all (stable) families in Table II will (slowly) converge, at very high 
orders, to approximately the same values, in a way again similar to what was observed in the simple geometric series 
example above discussed in Sec. IV. A and in the anharmonic oscillator case However, it appears that these PMS 
optimized solutions in Table II rather converge to a value of about ~ 2.2, i.e. close to but not equal to the exact 
large- iV result ci ~ 2.328. The reason for this slight discrepancy will be understood below. 

In contrast, the results in Table IV illustrate how good the IR approximation is, although from a perturbatively 
incquivalcnt approach, and indicates that the LDE series in Eq. H3.25|l do converge to the correct result. Pushing the 
LDE to higher orders, we obtain for Eq. H3.27|l . for instance at order 100, that most of the solutions are very close 
(within 0.1% error) to the exact result, the closest solution being 2.32834 ± 0.00132884/. 

In order to explain the convergence properties on more basic grounds, both for the exact 1/N or IR approximated 
series, one first observes that, similarly to the simple example studied in Sec IV. A above, all of the PMS optimization 
solutions corresponding to Tables II-IV, except for _F0, have I77I starting from relatively large values and then (rather 
slowly) decreasing as the LDE order k is increased, therefore reaching the boarder of the convergence radius of the 
original Su/r] series. More precisely, the convergence radius of the geometric series in Eq. (|3.27|l is immediately 
given by = l/(247r), so that, for u = 1, the original u/rj series (before LDE is performed) can only converge 
if \fj\ > l/(247r) ^ 0.013. [One expects similarly that the exact large- A'' original series in Eq. H3.21|l should have a 
convergence radius very close to this l/(247r) value, as we actually checked numerically by calculating e.g. = 
limi_^oo Ji+i/Ji to sufficiently high order i ^ 10"^]. On the other hand, the exact result should be recovered for 77 = 0, 
thus outside the convergence radius of the original perturbative expansion. Nevertheless, just as explained with the 
simple geometric series in Sec. IV. A above, this is compensated by the fact that the reorganized LDE series to order 
k, modify perturbative coefficients of order n with an extra damping factor of l/r[l + n/2] (see e.g. Eqs. (IB2p . (jB3p 
in Appendix B). As a result, this modifies the convergence radius so that smaller and smaller values of jryj can be 
reached, which basically explains e.g. the good convergence of the IR series to the right exact result shown in Table 
IV. Moreover, an important remark in view of the more interesting application to the finite N ~ 2 case, is to realize 
that most results produced by the families with real positive parts in Tables II-IV give reasonable values already at 
very low orders fc ~ 2,3, being for instance about 10% away from the exact result. We believe this is not at all 
a coincidence but simply reflects the faster convergence properties of the PMS solution within our prescription, as 
discussed in details with the simple example of Sec. IV. A, due to the presence of the linear tadpole term, which 
implied for the simpler geometric series case an exact cancellation for 77 — > of the first order term. For the actual 
IR geometric series the normalization of the tadpole term relative to the series is however completely different, but 
because \ff\ (247r)~"'^ as the LDE order p increases, one sees immediately that the first term of the expansion of 
the tadpole term for \fj\ (247r)-i is -N/{96tt'^), so that a similar cancellation still occurs, rendering again the 
convergence somewhat faster. 

Based on those general convergence properties considerations, we may also introduce a very simple criteria to select 
among the multiple PMS solutions as illustrated in Tables II-IV: we can consider the PMS solutions with the smallest 
1^1, yet such that |?7| is still within the convergence radius of the relevant series. Note that this criteria has also the 
advantage of automatically eliminating the negative FO solution, as it turns out that the latter always corresponds 
to the largest (and in fact, slowly increasing with k) \fi\ values, just like was the case for the simpler example of Sec. 
IV. A. Indeed, it is interesting to compare again the exact versus IR approximated 1/N series when replacing within 
their respective LDE expansion to order k this exact value of rj = R~^, instead of the PMS optimized fj values. This 
is illustrated in Table V, where it clearly appears that the IR approximated series behaves in a somewhat better way 
than the exact series. 

Now coming back to the exact large-A^ series case in trying to better understand the results in Tables II and 
III, we can alternatively study numerically the exact expression (|3.16|l directly, that is performing the integration 
numerically before expanding in u/r] series. In contrast with the latter expansion having (before LDE is performed) 
finite convergence radius as just discussed above, H3.16|l is of course defined for arbitrarily small 77. The behavior 
of this expression for u = 1 close to 77 = is shown in Fig. 2. As one can see, besides the exact result at 77 = 0, 
Eq. H3.1(j|l gets a real minimum ~ 2.06187 for very small 77 ~ 0.00235745. Now it appears clearly that all the LDE 
solutions in Tables II-III are disturbed by the presence of this real minimum very close to 77 = 0: namely in the 
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process of reaching smaller and smaller 77 values after the LDE-PMS procedure is applied, the PMS optimization 
solutions of the LDE at successive orders can be "trapped" into the first minimum reached. In contrast, notice on 
Fig. 2 that the IR approximated series has a similar real minimum but located exactly at r/ = 0, therefore it is not 
surprising that the PMS optimized LDE does converge correctly into this minimum in this case. This explains why 
the exact large- iV series converges to a value slightly different from the exact one, and is accordingly a weak point of 
the procedure. In view of this it is interesting to briefly compare the latter results with the alternative prescription 
such that the LDE-PMS procedure is applied on the same series, Eq. H3.21|l . but omitting the tadpole term linear in 

Tj: —Nr]* /{An) = {4>^)u'' — A(0^)c'^\ This corresponds to extremizing A(^c*^ directly (see Eq. ULyfl ). This is then 
similar to the recent LDE convergence studies of Braaten and Radescu j23l. Note that in this case, the additional 
extremum at very small 77 is removed. By selecting for simple illustration only the best PMS solutions (i.e. the ones 
whose real parts are the closest to the exact 1/A^ result ~ 2.328..), we obtain for ci the results shown in Table VI. 
As one can see, the results in this procedure do converge to the right result, however very slowly. On the other hand, 
one can also realize that the results at lowest orders fc ~ 3 — 4, are far away (more than a factor of 2) from the exact 
result. In fact, one should wait until about fc ~ 50 to have a reasonably good approximation. (Moreover, we note that 
may other solutions that are not shown in Table VI start to be very unstable at such high orders, some being several 
orders of magnitude away from the correct result). So, even if the series without tadpole term ultimately converges 
to the right result within this LDE prescription, it appears of not much practical use for the non-trivial N = 2 case, 
where we recall that only the first three perturbative coefhcients are known at the moment. This is to be contrasted 
with the results of Tables I-IV, where the fact that the lowest orders are already a good approximation is, as above 
explained, a consequence of the cancellation of the first order 0(77) term, and of the behavior of the PMS solutions 
as the LDE order k is increased. 

We will next discuss a method which fully exploits these scaling properties of the LDE-PMS solution at large orders, 
and which accordingly allows directly to resum the LDE perturbative series and to eventually further accelerate the 
LDE convergence, when applied to our BEG problem. 

C. Contour Integral Accelerated LDE Resummation Technique (CIRT) 



Having performed the usual LDE interpolation, with 77* = 771/1 — S and u Su, one obtains the physical quantity 
<i> expanded to order k. This procedure defines a partial sum that may be written as (see Eqs. (|3.13|l and (|3.21|l ) 

k 

$(^)(77*,fc) ^ Cn{u6r[f{r,*)r , (4.6) 

71=0 

where f{r]*) is a function of 77* whose form depends on the dimensionality of the physical quantity <i>. As we have 
already emphasized, the use of 77* is just an economical way to take into account the simple 5r]^ vertex. Usually, 
one expands 77* so that all the terms of order < k are present and the direct application of the PMS optimization or 
Fastest Apparent Convergence to this quantity, at 5 = 1, defines the standard LDE. On the other hand, one could 
be tempted to improve a perturbative series for which the highest order term contains {Su)'' and then expand 77* 
obtaining a higher order series. This may seem artificial at first because the contributions of order greater than k 
would come only from ^77^ insertions not taking into account, at the same order, new contributions which arise from 
the quartic Su vertex. However, by construction {4j| this procedure helps to accelerate convergence as it gives in a 
more direct and simple way the large LDE order behavior, as studied by the direct "brute force" LDE method in 
Appendix B. One may consider an all order resummation by observing that, for (5 — > 1, the partial LDE series is given 
formally, from the simple pole residues, as 

$ W (77, 5 -> 1) = ^ ^ dS $(77, u, 5) , (4.7) 

where the anticlockwise contour encircles the origin. Now, one performs a change of variables (see Ref. ^3 for the 
original application of this procedure) for the relevant (5 — > 1 limit, 

S=l-v/k, (4.8) 

together with a similarly order-dependent rescaling of the arbitrary mass parameter, r] ^ rj /c^/^, where the power 
1/2 is simply dictated by the form of the scalar mass interaction term, 77^ 0^ in Eq. 12. 6|) . This rescaling of the mass 
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parameter is of course consistent with what is obtained by a direct study of the large k behavior of the standard LDE, 
see Appendix B. 

For fc — > oo this resummation takes the final form of the replacement 

77* ^ Tyw^/^ , (4.9) 

followed by the contour integration 

$(fc-oo)(^ ^ 1) ^ J_ / ^ exp{v)'P{r,* ^ 77i>i/2) ^ (4.10) 



27ri / V 

where the "weight" exp(i;)/w originates from 

dS{l-6)-^ ^ -dv/v , (4.11) 

and 

lim (1 - v/k)-^-^ = exp(w) , (4.12) 

k — >oo 

while the original contour was deformed to encircle the branch cut Re{v) < 0. Here, one is initially dealing with a 
power series in {Su)[Su/ri*Y (cf. e.g. Eqs. H3.6|l . \S.H\\ ) and so, the use of 

dv exp (v) = 2TTi/r{-a) , (4.13) 

shows that the main effect of this resummation is to divide the original expansion coefficients at order-^*^*"*"^) by terms 
r(l + i/2) ^ (V2)! for large n. This damping of the perturbative coefficients at large order, as implied by this 
specific resummation, is fully consistent with what is obtained by a direct "brute force" resummation of the LDE 
series for large order k ^ 00, see e.g. Eqs. (jB2|) . (jB3|) in Appendix B. But such a damping is rather generic and 
was exploited recently in the completely different context of asymptotically free models 30] where it was shown to 
accelerate convergence of the LDE. When applied to the anharmonic oscillator, it is in fact (asymptotically) equivalent 
to the more direct LDE resummation with an order-dependent rescaling of the arbitrary mass, as employed in some of 
the Refs. |2^] to establish rigorous convergence of the LDE for the oscillator energy levels, that is itself an extension 
of the order-dependent mapping (ODM) resummation technique (see, Seznec and Zinn- Justin in Ref. 0|). In fact, 
this procedure can essentially suppress the factorial behavior at large orders of the perturbative coefficients generic 
in many theories, and convergence may be obtained even for series that are originally not Borel summable ,25^ .301. 
The above contour integral resummation is very convenient since it is algebraically simpler than the direct LDE 
summation (compare with Appendix B). In the present case, one expects fast convergence since the original series has 
no factorially divergent coefficients. 

Let us start by treating the standard closed form result of this contour integral accelerated resummation technique 
(CIRT) transforming Eq.ljS^B) into 



4^r(i/2) 3 ^r(i + V2)V J 

which clearly displays the coefficients damping. Then, by applying the PMS, Eq. 14.5|l . one obtains the results 
displayed in Table VII. This table shows that only the lowest order real parts of the positive families produce reasonable 
results but the values seem to deteriorate at high orders. Again, one observes that all positive families have similar 
behavior, starting from values close to ~ 2.000 the values decrease and then start to increase, as exemplified by Fl. 
Table VII seems to show that the disturbance due to the presence of the extra minimum of Eq. (|3.1t)|) : ~ 2.06187 for 
very small 77 ~ 0.00235745, reflected in the data of Table II, get amplified by the CIRT which resums the series. This 
effect will be further discussed in Sec. IV. E below. 

We shall now examine how the all order LDE summation can further improve the convergence of the series. To see 
that, we consider again the IR behavior (|3.25|l . that is sufficient to grasp the essential features by keeping all results 
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fully analytical. Applying thus the CIRT method to the simpler geometric IR series, Eqs. (|3.25|) and (|3.27|l . we obtain 
the result 
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y -1 
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Erfc{x) - 1 (4.15) 



with X = Nu/{A8ttt]) and Erfc{x) = 1 — Erf{x) is the standard error function (see also Appendix B). Next, upon 
using the well known asymptotic expansion of Erfc{x) for x ^ oo {i.e. 77 — > 0) |4l| (see Eq. (|B5p in Appendix B), 
we can again apply the PMS optimization procedure at given successive orders k. The PMS optimization generates 
the results shown in Table VIII. Note how the large-A^ result, ci = 87r/[3C(3/2)''/^] ~ 2.328 is exactly reproduced 
already at the lowest non-trivial order. Moreover, as one can see from that table, this real solution remains valid 
at any order in perturbation theory. Thus, the exact large-iV result is recovered from this exact CIRT resummation 
of the LDE, even if we have only used the IR approximation of the propagator. The convergence to the exact 1/A^ 
result in this alternative LDE implementation is extremely rapid upon using the CIRT. This also indicates that the 
convergence to the exact result can be independent of the details of the non-asymptotic behavior of the perturbative 
series coefficients, at least for the large- iV quantity here considered. 

An important feature of the CIRT-PMS results in Table VIII is that the asymptotic expansion in Eq. HB5() implies 
that we are now dealing with a series in rj/u, instead of the standard perturbative series in u/ry, which we started from, 
as in Eqs. (|3.13|) . (|3.21|) . or (|3.27|) . This is a consequence of the IR approximation, leading to the simple geometric 
series in u/f], Eq. (|3.27|) . whose all order CIRT form has an exact expression, Eq. (|4.15() . that can be re-expanded in 
a rj/u series. The reason why the convergence properties of such a series are much better than those of the original 
series should now be clear in view of the discussion in previous subsections: the original theory is to be recovered 
in the limit 77 ^ 0, which is clearly not in the (finite) convergence domain of the original u/rj series, while it is 
automatically inside the convergence domain of the all LDE order resummed series Eq. (|4.15|) . Moreover, another 
advantage of this reverted series approach is that we may bypass the need of PMS or other similar criteria: clearly 
the best approximation to the exact result will be simply given by the smallest rj values, irrespective of whether it is 
a solution of a PMS (or similar) criteria. 

D. LDE from Pade approximants 

For completeness we will also consider in this subsection the results obtained from a completely different resumma- 
tion of the relevant BEC series, based on the Pade approximant (PA) method We define, as is standard, a PA 
P[n, m] as a rational fraction of two polynomials of order n and m respectively in the relevant variable u/rj: 

where the perturbative coefficients and bj are obtained order by order by expanding Eq. (|4.m() up to order n + m 
and matching the resulting series with the original expansion. We recall that PA are generally useful when a series is 
known only up to the first few orders, as they can predict sometimes with a very good accuracy the unknown higher 
orders, and/or give very good resummation results. As far as the BEC series is concerned, our further motivation 
to consider PA is that it will allow to simply define an alternative (approximated) series in the inverse variable 77/ w, 
starting from the exact large N (or finite N as well) series Eq. (I3.13f) . As illustrated in previous subsection with 
the IR approximated series that has an obvious alternative expansion in rj/u, since the exact result is recovered for 
77 — > 0, we shall expect much better convergence properties from such an inverted series, as will be seen below. Note 
also that the PA method is largely independent of the previous methods, so that it can provide a further consistency 
cross-check of the numerical results. Another advantage is that the PA technique is immediately applicable as well 
to the finite N case, to be discussed below in Sec. V. A drawback of the PA method, however, is that they are 
not uniquely defined, since for a given perturbative expansion of order k one may consider a priory all possible PA 
with n + m = k. In order to thus limit somehow the number of possible PA without introducing much bias in our 
analysis, we will only consider resulting PA that can be expanded in powers of 77/u, for the reasons discussed above, 
which imposes that n < m. To illustrate in the most simple case the power of PA, let us first consider again the IR 
approximated geometric series Eq. (|3.27|) . but assuming that only the first order term of the u/rj series is known: 
{4>'^)iR ~ — (-^'7*)/(47r) + ((5wA^)/[(967r2)(87r)](— (5wiV)/(677*). We could then define an approximation of this series as 
follows 
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where the PA of order n + m — 1 is 



P[0, l]{u/ri) = [bo + b.SuN/ {67^)]-^ , 



(4.18) 



and a simple matching of the expansion of H4.17|l gives 6o = 1 and hi = l/(87r), such that the exact geometric series is 
in fact recovered, and can be of course expanded in the alternative form of a rj/u series. Though this example maybe 
too simple, we can expect that the more general PA that are constructed below to approximate the more complicated 
exact series, will have similarly good resummation properties. Typically at order 3 we have to consider P[l,2] and 
P[0, 3]. The PA results are shown in Table IX for the large N case. The order k designates in this case the order of the 
(re)-expansion in LDE power series of the PA, followed again by a standard PMS optimization. Higher orders k > 5 
are not shown but exhibit a very stable behavior with solutions almost similar to the lowest orders shown. As one can 
see, the exact result is often reproduced as a PMS solution, which is not so remarkable in the present 1/N case, as it 
simply means that rj = (for which value the PA is well-defined by construction) is a solution when applying the PMS 
to the PA. Clearly the best PMS solutions, if 77 0, are obtained from the smallest \ri\ ^ values. Though it appears 
clearly that the PA P[0, 3] give much better results than the P[l,2], note that the latter give solutions ^ 2.05 and 
^ 2.6 that are located not far below and above the exact value. The better P[0, 3] results maybe eventually explained 
by noting that the P[0, 3] should have much better resummation properties for 7/ — s- 0, having no numerator term in 
u/rj. 

Another advantage of the PA is that one may skip the PMS criteria and simply take the limit 77 ^ directly, 
in which case the exact 1/iV result is reproduced. Of course as concerns the large- case those results are only a 
consistency cross-check, since we started from the exact series anyway and thus only managed to define a PA such 
that its 77 limit is well-defined, in contrast with the naive u/rj series. But the very same PA procedure can be 
applied to the finite N case to be discussed in the next Section. 



One may now summarize the main features of this detailed large- A^ investigation so that the finite N application, 
in the next Section, can be carried out straight away. Regarding the family selection one notes, especially in Tables II 
and IV, that except for PO all families produce results that converge to approximate similar values at high orders, as 
in the anharmonic oscillator an geometric series cases. Therefore, choosing consistently PI in all the BEG applications 
appears to be an appropriate choice since this family is also the only one which allows predictions at any order {e.g. 
in a computation involving only low orders in S). Again, this is consistent with the observations drawn from the 
geometric series analysis. 

To understand completely the reason for convergence it is particularly convenient to examine the formal expressions 
of the large LDE order behavior, i.e. for k 00, as derived in details in Appendix B. What happens, as already 
explained in Sec. IV. B, is that when the LDE order is increased, the numerical PMS optimal solutions \fj\ tend to be 
smaller and smaller and thus to reach the boarder of the convergence radius of the original perturbative series. But 
at the same time, from the study of the asymptotic behavior of the series at large LDE orders we see that the main 
effect of the rearranged LDE series is to provide an extra damping factor l/r[l + 7t/2] as well as a scaling factor A:"/^ 
in its order n coefficient. So, after redefining conveniently r] — > f/k^^'^ , the new optimized 77 values can tend to zero, for 
which the new LDE-resummed series is now more and more convergent (see 6.17. Eq. H4.15|l ) and leads to the correct 
77 ~ result. However, in the process of reaching smaller and smaller values, it may happen that another non-trivial 
extremum, if present, is first met, in which case the convergence is disturbed or slowed down by the presence of this 
other extremum and the 77 — may no longer be reached. This is precisely what happens in the case of the exact 1/A^ 
series, as illustrated above in Tables II and III, with this problem worsening as the series is resummed, as shown by 
the results in Table VII, for the GIRT resummation applied to the exact large- A^ perturbative series. In contrast, by 
using the IR series approximation the disturbing extremum that now approaches 77 7^ is removed and this problem 
is not present there. Finally, we mention that the recent LDE convergence studies of Braaten and Radescu zj\ follow 
the same lines of the standard LDE calculations performed by us in Sections III.B and IV. B. The major difference 
being that those authors prefer to extremize lS.{4P')'"p directly (see Eq. ULyfl ). In practice, the difference in between 
those physical quantities amounts to the tadpole term —Nrj* /{An) — {(jP')u^ — A{(jP)i^\ that appears in Eqs. (|3.13|) 
and (|3.21l) . In the case of the standard LDE-PMS application, like the one shown in Table II, their numerical results 



E. Brief summary and discussion of the convergence properties 
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are similar to ours at very large orders, exhibiting ultimately convergence, but do not give a good approximation at 
the lowest orders. The reason for this faster convergence when the tadpole term is present is the cancellation of the 
leading terms for ^ as explained in detail above. Our simple geometric series investigation has also shown the 
crucial role, regarding convergence, played by this type of linear term. Moreover, these type of loop terms are also at 
the origin of good convergence properties observed in many other applications (2^ i3^J . 



V. THE FINITE TV CASE WITH RESUMMATION OF THE LDE 



Let us now turn our attention to the finite N results and, especially, to the improvement of the LDE within this 
limit by using the GIRT resummation method discussed in the previous Section. The insight gained in the detailed 
large- TV study will prove to be very useful to understand the structure of optimized results and to select the appropriate 
solutions. For arbitrary N, the quantities {4>^)''u^ as well as 5r^^^ have been evaluated in details, up to order-J'', in Ref. 

|18| . The contributing diagrams evaluated in that reference for the 5 perturbative expansion of {4>^)'^u'^ to order-5* are 
shown in Fig. 3. In Ref. [ig all the integrals appearing in those diagrams were obtained with the type of perturbative 
calculation discussed in Subsec. III.A with the multidimensional Feynman integrals calculated with VEGAS. Those 
terms were also later been obtained by Braaten and Radescu in a different way, by reducing the multidimensional 
integrals to one-dimensional ones in some of the terms of fourth order in 5 and, more recently, Kastening in Ref. 0| 
has also revised those numerical results, obtaining more precise numerical results for the integrals. From the results 
originally obtained in Ref. 18] and using the corresponding updated, higher precision coefficients evaluated in |44|. 
the terms contributing to to order-J**, shown in Fig. 3, are given by 



(^*)2 (47^)6 108 ^ ■ J ' ryS (i8)2(47r)7 



where, as before, 77* ~ 77%/! — 5 must be expanded accordingly to take into account all contributions shown in Fig. 3. 
The symmetry factors appearing in Eq. can be found e.g. in ]^^- The symmetry coefficients clearly show that 

this perturbative expansion is valid for any N which means that, up to order-(5^, the large- A'^ results Eqs. H3.13|l and 
(|3.21|) may be recovered, within numerical error bars of about 2%. This evaluation is easily done by considering u to 
be of order 1/iV so that only {uN)'' terms are retained together with the first, u-independent term. Fig. 3 illustrates 
well how the LDE mixes, at a given order, diagrams which normally appear at different orders in the 1/TV expansion. 
By setting N = 2 one then gets the more compact form 



mi'^--'^+Suj2K,^~^:^y + OiS') , (5.2) 

where the coefficients are Ki = 3.22158 x 10"^, K2 = 1.51792 x 10^*^, K3 = 9.66514 x 10"®. It is worth remarking 
that the coefficients Ki and K2 obtained from the numerical results of Ref. agree with these results, also obtained 
later by the authors of Refs. 23, 44]. At the same time the coefficient used here, which was obtained in Refs. 
[23,01, differs by about 10% from the one that would come from the results of Ref. 0. We in principle can trace 
this difference by the fact that five non-trivial graphs with five loops enter the evaluation of (see Fig. 3), which 
require some careful calculation^. 

Then, by applying the standard PMS, Eq. 12.21) . to Eq. I|5.2f) one obtains three families, shown in Table X, in 
agreement with Ref. |0]. Turning to the GIRT resummation of Eq. 15.2|l one proceeds as in the large- case (see Sec. 



^ We thank B. Kastening for pointing out to us the correct values for the five-loop diagrams and for discussions concerning the difficulties 
in evaluating them. 
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IV. C). By applying the CIRT improved in the context of the PMS optimization to Eq. H5.2|l one obtains three other 
famihes, also shown in Table X. In fact, the previous large- analysis strongly suggest that also here the first family 
with positive real parts should be the relevant family concerning our N = 2 predictions. Indeed, if we assume that 
the large order behavior of the actual N = 2 series coefficients should not be drastically different from the analogous 
large- series, all the fast convergence and scaling properties that were discussed in Sec. IV should be approximately 
valid also for N = 2. Then, the first negative family in Table X is easily eliminated by the same criteria as in the 
large- case, because it again always corresponds to the largest |^| and does not exhibit any trend towards smaller \fj\ 
values as the LDE order k is increased. Similarly, we also notice that the Fi family in Table X has Re(77) substantially 
smaller than F2, while we expect the exact result to be for ry ^ 0. Moreover, due again to the presence of the tadpole 
term in our procedure, from the analysis in Sec. IV we can expect that our results, though intrinsically limited for 
iV = 2 to the first four LDE orders, should be nevertheless already a reasonably good approximation. 

A further crosscheck of the consistency of our results without any knowledge of the exact higher order coefficients for 
the case A'' = 2, is the stability of the result when replacing these unknown perturbative coefficients by a well-defined 
approximation. This is the result exhibited in Table XI, where the unknown order Ki with i > 4 were replaced by 
the corresponding coefficients of the IR approximated series, Eq. H3.27|l . As one can see, the stability of those results 
is quite remarkable, and one can even observe the slow convergence of the standard LDE-PMS to the CIRT result. 

The physically meaningful real part of our order-J'* CIRT improved result ci ~ 1.19 can then be compared with 
the recent Monte Carlo estimates ci = 1.32 ± 0.02 and ci = 1.29 ± 0.05. Note that the standard order-5'' PMS result, 
Ci ~ 1.53 is also a satisfactory estimate. We note that the CIRT and ordinary PMS results just bound the Monte 
Carlo estimates from below and above, respectively. 

Finally, in Table XII we show the results for ci obtained by Pade approximants (PA), as discussed in the last 
Section. As expected, at first non trivial LDE order (order 3), only the 1/A^ solution is found from the PA, because 
only 77 = is a PMS optimization solution. As one can see, the non-trivial results at higher orders (5-10) of the LDE 
expansion in Table XII are nicely consistent with what is independently obtained from the standard LDE and CIRT 
results shown in Table X. Note however that the PA P[0,3] are in better consistency than the P[l,2], which in fact 
only give result very similar to the ones in the large- A?^ case shown in Table IX. 

The contributions to Sri^'' = — l]ren(0), which enter in the derivation of the constant Cj, Eq. (|1.6|l . have also been 
explicitly evaluated in Ref. 18] and again we refer the interested reader to that reference for the details and we show 

below only the final result for the renormalized, scale {M) dependent Sri'^^ obtained in that reference: 



Btt V 3 J (47r)2 18 



In ( — 1 - 0.59775 

1]* 



,ou3(A^ + 2)2^ , „u3 (16 + 10Ar + Ar2) 

r;2 6(47r)6 108 ^ ^ 



,^(40 + 32A. + 8^2 + ^3) _ ^(44 + 32iV + 5A.2) 

r;2 (47r)6648 ^ ^ 772 (47r)6324 



which, for A^ = 2, becomes 
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^'-;A, + 5^^^A, + 0{5') , (5.4) 



where A2 = 1.40724 x lO'^, A3. = 8.50859 x 10"^ and A^ = 3.52299 x 10"*^. The application of the PMS optimization 
to Eq. (|5.4|l reproduces the same results as obtained in Ref. |0 which are ^ 101.4, ^ 98.2 and ^ 82.9 from second 
to fourth order respectively. 

At the same time, treating Sr'^f^ with the CIRT (from Sec. IV. C) one obtains the result^ Re [r'^f^ [M — u/3)] = 
0.0010034^2 which, together with the CIRT improved {4P')u'' result, and Eq. ifL^ . leads to (with errors estimated 



^ Note that the scale M = u/3 was originaUy chosen in the Monte Carlo applications 
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from the integrations performed in 1^ with Vegas) = 84.9 ± 0.8 whereas the Monte Carfo result is = 75.7 ± 0.4 
i^J. To our knowledge, these are the only analytical predictions for this coefficient to the present date. For consistency, 

note that the optimization of (5rc^^ , including the selection of solutions, has also been performed according to what 
was done for (0^)1''^ 

Note that we have not attempted to examine the infrared behavior of the finite N case since its series is much more 
complicate than the large- one. However, our previous large- investigation shows that the LDE works well, even 
for the standard series, already when considering only the lowest order terms. Here, only these lowest order terms 
were computed so that our results, up to order-J^, can be considered good estimates even if one knows that the whole 
procedure may get spoiled at very high orders. 



VI. CONCLUSIONS 



We have investigated how the LDE followed by a standard PMS optimization performs in the non-trivial case of 
phase transitions of interacting homogeneous dilute Bose gases described by an effective three dimensional (f)'^ field 
theory. This nonperturbative method has been recently employed in Refs. [l3 ^-^d in Ref. 23] to determine the 
critical temperature for such a system, giving good numerical results. One advantage is that the formal calculations 
are performed exactly as in the perturbative case. This means that, at each order, one deals with a very reduced 
number of contributions which are not selected according to their topology (like the number and type of loops). 
Therefore, the method is valid for any finite value of N. To handle ultraviolet divergences, the renormalization 
program is implemented as in the usual perturbative way. Also, an arbitrary mass parameter consistently introduced 
by the method avoids any potential infrared problems. 

The convergence properties, including rigorous proofs, of those nonperturbative methods have been studied in 
quantum mechanics |25l l26l l27l |2^ and more recently in quantum field theory |30l | . However, despite the many 
successes obtained with the LDE in different applications, the convergence study in the BEG case poses new challenges. 
One of the reasons is that it is difficult to establish simple analytical links between the LDE and other nonperturbative 
methods at the one-loop level, since those terms do not contribute at the transition point. This could rise the suspicion 
that the results obtained in Refs. 17, 18], for the realistic N = 2 case, are just a numerical coincidence. The exact 
value for the linear coefficient ci which appears in the critical temperature of interacting homogeneous dilute Bose 
gases is still unknown |^ although much progress has been recently made concerning its determination 13]. Here, 
our aim was to prove the reliability of the recent LDE results for ci and .18] through a detailed analysis of the 
LDE convergence properties. 

We have started our convergence study by considering the effective BEC model in the large- where the results 
obtained by Baym, Blaizot and Zinn- Justin (ci ~ 2.328) [i3 can be considered "exact". We have performed the usual 
perturbative LDE evaluation of {(jp)u'^ in different ways which allow for numerical accuracy checks. By considering 
the asymptotic infrared behavior of the propagator we have obtained a simpler series with exact coefficients which 
allows for a fully analytical analysis. Before tackling the optimization of the BEC series we have investigated how 
the procedure works considering a simple geometric series. The insight gained during this exercise proved to be 
very important in understanding the family structure of optimal solutions with regard to convergence properties. 
Then, the standard series and its infrared limit have been optimized with the PMS criterion, Eq. (|2.2|l . leading 
to reasonable results in the standard case (see Table H). At the same time, the numerical results produced by the 
same optimization applied to the infrared series display better convergence properties, as shown in Table IV. This 
also shows that the optimization procedure is rather sensitive to the actual form of the perturbative expansion series 
coefficients, in particular at low perturbative orders. Indeed, we emphasize again that including the linear tadpole 
term in the LDE-PMS procedure is crucial for a faster convergence and to obtain an already very good approximation 
at low perturbative orders, though both procedures tend to similar results at very high orders. 

We have then presented an efficient all order resummation technique, similar to the one used to prove the LDE 
convergence within quantum field theories at zero temperature 30] . This LDE resummation method takes advantage 
of contour integration techniques which allows to resum more directly the series and thus to accelerate convergence. 
Applying this contour integral resummation technique (GIRT) to the exact large- TV perturbative LDE series seems 
however at ffist to amplify the numerical instabilities generated through numerical optimization. This problem 
becomes more severe as one moves to higher perturbative orders as shown in Table VII. On the other hand, applying 
the GIRT to the infrared series exactly reproduces the large- iV value ci = 2.328 already at the first non-trivial order. 
This solution has no complex parts and remains valid as one goes to higher orders while all other families of solutions 
display good convergence properties having real parts that are numerically very close to the exact value (see Table 
VIII). In summary, this extensive analysis has shown that, for this type of series, the optimization procedure and 
convergence rate may be influenced by the actual values of the first few perturbative order coefficients. We have 
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shown that, for this effective BEC model, the simple infrared series retains all the nonperturbative information and 
that the LDE, augmented with the CIRT, performs rather well. Finally, we have seen that all families, display a very 
similar structure and will predict approximately the same values at high perturbative orders. We also have shown 
how the family which produces negative ci values is easily eliminated because it does not correspond to the expected 
trend towards small \ri\ values of the PMS solution. The very same criteria allows to single out the values generated 
by the first family of positive solutions, both for the large- iV or finite N case. Finally, the results obtained by a 
different resummation based on Fade approximants, aimed as an alternative to define the relevant 77 —s- limit from 
the perturbative series in u/rj prior to LDE, appear quite consistent with the latter ones. 

More formally, our present study investigated in some detail the large order behavior of the LDE (e.g. in Sec. IV. C 
with the CIRT method or alternatively in Appendix B) , from which we can also point out interesting analogies as well 
as differences with the LDE convergence properties in quantum mechanics for the anharmonic oscillator 25, 26, 27, 2^. 
The latter is described by a (one-dimensional) scalar theory with a cj)'^ interaction term, and as is well-known its energy 
levels have perturbative expansion coefficients that are factorially growing at large ordersls^. Nevertheless, as already 
mentioned in Sec. II, the LDE can converge essentially because the PMS optimized solutions behave like a rescaling of 
the mass parameter with perturbative order which can compensate the factorial behavior at large orders. In contrast, 
the relevant BEC perturbative series here considered in Sec. Ill have a finite convergence radius, such that no explicit 
rescaling of the mass parameter should be necessary in principle for convergence. Nevertheless, what the LDE followed 
by PMS optimization is performing is to enlarge the original series convergence radius, thus qualitatively similarly in 
this respect with the oscillator. These properties of the LDE are best exploited by the CIRT more direct resummation 
method. 

The final part of the work was devoted to the realistic finite N case for which only the standard series with coefficients 
numerically obtained is available. In practice, here only the first low orders contributions could be evaluated and the 
comparisons performed in the large- TV show that VEGAS produces, in this case, accurate coefficients which should 
not completely spoil the optimization procedure. For consistency with the large- A'^ case we have considered only the 
real parts of the first family of positive solutions as the relevant ones. Applying the CIRT to this case has improved 
the recent order-J^ results of Ref. jl^, generating ci ~ 1.19 and Cp ^ 84.9, which are about 9% smaller and 11% 
higher, respectively, than the recent lattice Monte Carlo estimates '^^^^ cannot expect to make 

a definitive analytical prediction, for those coefficients, from a calculation involving only a handful of contributions. 
Nevertheless, the agreement between our improved analytical LDE results and the recent numerical Monte Carlo 
results is quite impressive. Moreover, the consistency of our N = 2 results as obtained from a different resummation 
method (Fade approximants) used prior to the LDE procedure is also noticeable. Our results seem to support the 
fact that, analytically, the leading contributions to ATc for the BEC case studied here can be obtained by resumming 
typical leading and next to leading 1 /N type of graphs as in Refs. [13 and [13. The better LDE numerical values 
may be due to the mixing of such contributions since in reality, N — 2. 

In summary, our detailed convergence study together with the improved optimization procedure results show the 
potential of the LDE to tackle nonperturbative calculations in field theory at critical points. We have explicitly shown 
how meaningful non perturbative results for the BEC problem can be obtained in a consistent fashion which also 
works for a general case such as the simple geometric series analyzed in the paper. 
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APPENDIX A: REVIEW OF THE ORIGINAL LARGE-A CALCULATION. 

Let us briefiy recall how the exact large- derivation of ci was performed in the original calculation, Ref. |lClj |. in 
order to exhibit the differences with the LDE (exact or approximated) evaluation as performed in Sec. III. Conside ring 
the original theory described by Ea. (|1.8|) at the critical point, the Hugenholtz-Pines (HP) theorem imposes fA Il4j 
Tc = — S](0), where S(0) represents the field self-energy with zero external momentum. Then, one has a massless 
propagator and the appearance of IR divergences has to be carefully dealt with. In Ref. T^, after applying the HP 
theorem, the relevant expression for A{(jP) reads 



(27r)3 



1 



p2 + E(p) - S(0) p2 



(Al) 



27 



where l/p^ represents the term with no interaction, {(fp')u^o, to be subtracted according to the discussion in Sec. II 



(compare e.g. with Eq. 1)1. 9|) ). and 
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with the "dressed" (resummed) scalar propagator (see Fig. 1) 
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represents the basic one-loop (massless) bubble integral. At the relevant (next-to-leading) 1/N order one obtains the 
basic expression to be evaluated as 
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where the subtraction term cancels out with the first order term belonging to the 1/A^ expansion of the term 
[p^ -(- S(p) — I](0)]~^, in Eq. (|A1|) . Remark that, after applying the HP theorem, all loop integrals in l|A5|l involve 
massless scalar propagators. Accordingly, both integrals in Eq. (|A5|) are IR divergent at intermediate steps, although 
the final physical result A(0^) should be an IR (and UV) convergent quantity. Also, the integral over k is not 
(absolutely) UV convergent: it has superficially a logarithmic UV divergence. Thus, as emphasized in Refs. [lol[Tl| . 
one should be very careful to correctly regularize these integrals before doing standard manipulations, like typically 
exchanging the order of the two integrations in Eq. (|A5p . The authors of Ref. ^3 chose to work in dimensional 
regularization, which takes care of both, UV and IR divergences. We thus now summarize the main steps of this 
calculation. First, one integrates over p, so that the second term of the bracket in Eq. l|A5p vanishes, since / d'^p/p'^ — 
in dimensional regularization. The integration of the non- vanishing first term in the bracket of Eq. (|A5p gives a result 
whose behavior for d ^ 3 is essentially given by ^ l/[T{d — 3)] , that would, naively, give zero as a result. However, 
it combines with a simple pole in d — 3 given by the next k integral (see below) . More precisely, integration of Eq. 
[51 gives (omitting factors that are regular for d 3): 



d'p 1 _ 1 r(d/2-l) TT 



{2^^ p-^ {k + pY (47r)'i/2 r(d-3) sin(7rd/2) 



(A6) 



where the space-time dimension d is kept arbitrary, for the moment. One has next to deal with an integral over k of 
the generic form (omitting again nonessential constant overall factors) 

{2^r ^+B{d)k'i-^ ' ^ ^ 

where B{d = 3) = 1/8. Next, to evaluate Eq. I|A7|I . we make the following change of variable: = k~'^ where 
e = 4 — d, which gives 

ddk^^p-^^^+y^WP . (A8) 
Then using the basic dimensional regularization formula: 

(27r)'i (p2 + i?2)a - (2^)d/2r(d/2)r(a) '^/^J J ' l^yj 
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we obtain after some algebra 

„2(d-3)/(4-d) ( .4 , 

't'- (47r)'*/2 (4-d) r(d/2)sin(27r(3-d)/(4-d)) \NuJ ' ^ ' 

where, for d — > 3, sin' "^[27r(3-d)/(4-d)] combines with r^i(d-3) to give a finite result^ proportional to u. Putting 
back all overall factors, one finds 

= , (All) 

which is simply related to ci via Eq. Ijf .5|l . in agreement with Rcf. [Tfil |. 

APPENDIX B: LARGE ORDER BEHAVIOR OF STANDARD LDE AND PMS 

In this appendix we briefly analyze the large order behavior of the standard LDE, in order to exhibit some generic 
properties of the LDE-PMS optimization solutions, as well as the link with the more direct CIRT method considered 
in Sec. IV.B. 

Considering either Eq. (|3.2HI or H3.27|l with ry* = 77(1 — (5)^/^, the result of its expansion to order k in S followed 
by 5 — > 1 can be written formally as (the coefficients Kn here refer indistinguishably to either the J„ of Eq. H3.21|l or 
the G^ of Eq. (|T77|) ) 

r[3/2] uN'^^ ( uNY"-^'^ r[l-,^/2] 



^~ 4. t'o -!r[3/2-n]+ 3 ^^^^^ 6 J ^/ q\T[l-n/2-q] ^^^^ 

where the last sum originates from the expansion of (1 — J)""/^ and the upper limit of this sum takes into account 
that at order k of the S expansion, there is a term coming from u{u/'q)"- ^"+^ u{u/r])"- and a term from (1 — (5)""/^. 
The sums can be performed analytically to give 

A/^2N(fc), ^ Nrj ^ (-1)'= "^^.^ f r[-n/2 + fc] 

^^'^^^ - - 4^^ r[i/2 - fc] r[i + fc] + - g ^" l"^ j r[-nU]r[iW2] ^^^^ 

The expression (jB2p as it stands is particularly convenient to be optimized with respect to rj at arbitrary order k, 
leading to the results shown e.g. in Tables II- VI, X, XI, depending whether one takes the exact 1/N, IR approximated, 
or exact N — 2 values of the relevant perturbative coefficients Kn, given in Eq. 1)5.2(1 . 

The large LDE order behavior, for k 00, of expression (|B2p can also be analyzed, to give: 



fc^oo 47r V7FfcV2 3 y y r[l + n/2] 

where we used, to obtain Eq. l|B3p . the well-known properties of the Gamma functions, such as r[2z] = 
22^-i/2r[z]r[l/2 + z]/y(27r), r[z]r[l - z\ sin(7rz) = tt and the Stirling asymptotic behavior V[b + az] ~ 

Expression l|B3p clearly suggests to rescale for convenience the arbitrary mass parameter according to 77 — > fjk^^'^. 
Of course, after such a rescaling the relevant limit is again 77 ^ 0. After such a rescaling, one obtains 

A(02)M(r7,.)^-^^ + ^y^^^ f^^X (B4) 



47r r[I/2] 3 ^^r[I + n/2] V 677 J 



Note a misprint in Eq. (27) of ref. llOl . where the relevant term sin[27r(3 — d)/{4 — d)] reads sin[n{d — 2)/{4 — d)]. 
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which, as expected, agrees with the CIRT "direct" LDE resummation result Eq. (|4.14|l . 

More precisely, for the simpler geometric series cases Eq. (14.21) or Eq. H3.27|l . corresponding thus to (up to an overall 
factor) Kn = 1, the sum in Eq. (|B4p can be further performed exactly, to give: exp(a;^) Erfc[x) — 1 with x = u/rj 
(or X = uN/{A%TTrj)) for Eq. (|4.2|l (respecti vely for Eq. H3.27II '). Finally the asymptotic expansion of Erfc{x) for the 
relevant limit x — *■ oo (equivalently 77 — > 0) |4l| : 




was used at different stages in Sec. IV, for instance to examine the behavior of the LDE series when the linear tadpole 
term is included in the procedure. 
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= xx>-<x 



FIG. 1: The Feynman graph for the relevant quantity {4>'^)u at 1/N order, with the resummed propagator (dotted lines). 
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FIG. 2: Comparison between the naive IR, UV propagator approximation Eq. H3.25|l . 1)3.26^ . and the exact numerical 
integration of Eq. 13.161 . before subtracting the spurious contribution A^ii/(967r^). 
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FIG. 3: All diagrams contributing to {4>^)'i^'' at the critical point. The black dot represents the Srj'^ insertions. 



TABLE L LDE of the simple alternated geometric series, Eg. (14. It . PMS optimization results, at different LDE orders k. 



K 




CI 


CO 






JO 


1 




-^2 










2 


1.69743 


-1.12996 














±0.167336/ 










3 


1.81088 


-1.00296 


-1.11978 












±0.129109/ 










4 


1.87082 


-0.959082 


-1.06361 












±0.076116/ 


±0.0465028/ 








5 


1.90768 


-0.950095 


-1.02417 


-1.05508 




— 






±0.03648/ 


±0.0499/ 








6 


1.93258 


-0.954501 


-1.00183 


-1.03612 










±0.01086/ 


±0.0407/ 


±0.019/ 






7 


1.95052 


-0.9630 


-0.99069 


-1.01992 


-1.03154 








±0.00447/ 


±0.02932/ 


±0.02416/ 









i.yD4Uo 


-u.y / iy 


-U.yoDi 


-i.UUoO 


-i.Uzoio 








±0.01304/ 


±0.01929/ 


±0.023/ 


±0.0096/ 




9 


1.97462 


-0.9797 


-0.9851 


-1.0012 


-1.01516 


-1.0205 






±0.0174/ 


±0.0114/ 


±0.0195/ 


±0.01346/ 




10 


1.9831 


-0.9861 


-0.9859 


-0.9968 


-1.0088 


-1.016 






±0.019/ 


±0.0056/ 


±0.0154/ 


±0.0141/ 


±0.0055 



TABLE II: Standard LDE at large- A*', Eq. 13.2111 . PMS results for ci, at different orders k, obtained with all families. 



K 


r U 


171 
r L 


ZT'O 

r z 


t 6 


TP A 

r 4 


r 


t D 


t I 


r 


Z7n 

t y 


r iU 


2 


-2.163 


2.163 




















3 


-2.698 


1.879 
























±0.169/ 




















4 


-2.945 


1.713 


1.962 






















±0.061/ 




















5 


-3.087 


1.642 


1.913 






















±0.077/ 


±0.029/ 


















6 


-3.179 


1.620 


1.870 


1.935 




















±0.196/ 


±0.008/ 


~ 
















7 


-3.244 


1.622 


1.843 


1.928 




















±0.293/ 


±0.067/ 


±0.001/ 
















8 


-3.292 


1.636 


1.829 


1.921 


1.935 


















±0.369/ 


±0.128/ 


±0.025/ 
















9 


-3.329 


1.654 


1.825 


1.913 


1.936 


1.946 
















±0.429/ 


±0.184/ 


±0.062/ 
















10 


-3.358 


1.675 


1.828 


1.908 


1.938 


1.950 
















±0.476/ 


±0.233/ 


±0.102/ 




±0.027/ 












11 


-3.382 


1.696 


1.836 


1.907 


1.940 


1.950 


1.966 














±0.515/ 


±0.276/ 


±0.139/ 




0.056/ 












12 


-3.402 


1.715 


1.845 


1.909 


1.941 


1.950 


1.972 














±0.546/ 


±0.313/ 


±0.174/ 




±0.086/ 


±0.025/ 










13 


-3.419 


1.733 


1.856 


1.913 


1.942 


1.950 


1.974 


1.983 




__ 








±0.571/ 


±0.344/ 


±0.206/ 




±0.114/ 


0.050/ 










14 


-3.433 


1.750 


1.868 


1.919 


1.942 


1.952 


1.976 


1.988 












±0.592/ 


±0.372/ 


±0.234/ 




±0.141/ 


±0.075/ 


±0.023/ 








15 


-3.445 


1.765 


1.880 


1.926 


1.943 


1.956 


1.978 


1.992 


1.997 










±0.610/ 


±0.395/ 


±0.260/ 




±0.166/ 


±0.098/ 


±0.045/ 








16 


-3.456 


1.780 


1.891 


1.934 


1.944 


1.960 


1.981 


1.995 


2.002 










±0.625/ 


±0.416/ 


±0.282/ 




±0.189/ 


±0.119/ 


±0.066/ 


±0.021/ 






17 


-3.466 


1.792 


1.903 


1.942 


1.944 


1.966 


1.984 


1.998 


2.006 


2.010 


— 






±0.638/ 


±0.433/ 


±0.302/ 




±0.210/ 


±0.140/ 


±0.085/ 


±0.040/ 








Q /I 7/1 


1 on /I 
i.oU4 


i.yio 


i.you 


1.940 


i.y / z 


i.yoo 


Z.UUU 


Z.UU9 


z.Ui4 








±0.649/ 


±0.449/ 


±0.321/ 




±0.229/ 


±0.159/ 


±0.104/ 


±0.059/ 


±0.019/ 




19 


-3.482 


1.815 


1.924 


1.959 


1.945 


1.978 


1.992 


2.004 


2.012 


2.018 


2.020 






±0.658/ 


±0.463/ 


±0.337/ 




±0.246/ 


±0.177/ 


±0.121/ 


±0.076/ 


±0.036/ 




20 


-3.488 


1.825 


1.933 


1.967 


1.945 


1.984 


1.997 


2.007 


2.016 


2.021 


2.025 






±0.666/ 


±0.475/ 


±0.351/ 




±0.262/ 


±0.193/ 


±0.138/ 


±0.092/ 


±0.053/ 


±0.017/ 



TABLE III: Standard LDE at large- iV, Eq.lEHJ. Three families of PMS and Fastest Apparent Convergence (FAC) results for 

Cl. 



k 


Fl - PMS 


Fl - FAC 


F2 - PMS 


F2 - FAC 


F3 - PMS 


F3 - FAC 


2 


2.163 


2.498 










3 


1.879 ± 0.169J 


1.884 ± 0.274/ 










4 


1.713 ± 0.061J 


1.660 ± 0.097/ 


1.962 


1.996 






5 


1.642 ± 0.077J 


1.593 ± 0.078/ 


1.913 ± 0.029/ 


1.917 ±0.0462/ 






6 


1.620 ± 0.196/ 


1.585 ± 0.213/ 


1.870 ± 0.008/ 


1.858 ± 0.003/ 


1.935 


1.941 


7 


1.622 ± 0.293/ 


1.600 ± 0.314/ 


1.843 ± 0.067/ 


1.826 ± 0.063/ 


1.928 ± 0.001/ 


1.928 ± 0.005/ 


8 


1.636 ± 0.369/ 


1.622 ± 0.390/ 


1.829 ± 0.128/ 


1.813 ± 0.130/ 


1.921 ± 0.025/ 


1.915 ± 0.0213/ 


9 


1.654 ± 0.429/ 


1.647 ±0.448/ 


1.825 ± 0.184/ 


1.811 ± 0.189/ 


1.913 ± 0.062/ 


1.905 ± 0.061/ 


10 


1.675 ± 0.476/ 


1.671 ± 0.493/ 


1.828 ± 0.233/ 


1.817 ± 0.241/ 


1.908 ± 0.102/ 


1.899 ± 0.102/ 


11 


1.696 ± 0.514/ 


1.694 ± 0.529/ 


1.836 ± 0.276/ 


1.827 ± 0.284/ 


1.907 ± 0.139/ 


1.898 ± 0.141/ 


12 


1.715 ± 0.545/ 


1.715 ± 0.558/ 


1.845 ± 0.313/ 


1.839 ± 0.321/ 


1.908 ± 0.174/ 


1.901 ± 0.178/ 


13 


1.733 ± 0.571/ 


1.734 ± 0.582/ 


1.856 ± 0.344/ 


1.852 ± 0.353/ 


1.913 ± 0.206/ 


1.906 ± 0.210/ 


14 


1.750 ± 0.592/ 


1.752 ± 0.601/ 


1.868 ±0.371/ 


1.865 ± 0.380/ 


1.919 ±0.234/ 


1.913 ± 0.239/ 


15 


1.765 ±0.610/ 


1.767 ±0.618/ 


1.880 ±0.395/ 


1.877 ± 0.403/ 


1.926 ± 0.260/ 


1.921 ±0.265/ 


16 


1.780 ± 0.625/ 


1.782 ± 0.631/ 


1.891 ±0.416/ 


1.890 ± 0.423/ 


1.934 ±0.282/ 


1.930 ± 0.287/ 


17 


1.792 ± 0.638/ 


1.795 ± 0.643/ 


1.903 ± 0.433/ 


1.902 ± 0.440/ 


1.942 ± 0.302/ 


1.939 ± 0.308/ 


18 


1.804 ± 0.649/ 


1.807 ± 0.654/ 


1.913 ± 0.449/ 


1.913 ±0.455/ 


1.950 ±0.321/ 


1.947 ±0.326/ 


19 


1.815 ±0.658/ 


1.817 ±0.663/ 


1.924 ± 0.463/ 


1.923 ± 0.468/ 


1.959 ± 0.337/ 


1.957 ±0.342/ 


20 


1.825 ± 0.666/ 


1.827 ± 0.670/ 


1.933 ± 0.475/ 


1.934 ± 0.480/ 


1.967 ±0.351/ 


1.966 ± 0.356/ 



TABLE IV: Infrared LDE at large- iV, Ea. iTTfl . PMS results for ci, at different orders k. 



K 


-T U 


171 
r L 


ZT'O 

r Z 


r 6 


TP A 

r 4 


t 


r D 


r 1 


t 


r 9 


171 n 
r iU 


2 


-2.852 


2.852 




















3 


-3.577 


2.444 
























±0.276/ 




















4 


-3.910 


2.244 


2.482 






















±0.200/ 




















5 


-4.100 


2.184 


2.397 






















±0.097/ 


±0.079/ 


















6 


-4.223 


2.184 


2.333 


2.397 




















±0.020/ 


±0.081/ 


~ 
















7 


-4.309 


2.205 


2.298 


2.369 




















±0.028/ 


±0.060/ 


±0.032/ 
















8 


-4.372 


2.232 


2.283 


2.342 


2.366 


















±0.055/ 


±0.037/ 


0.040/ 
















9 


-4.214 


2.256 


2.279 


2.324 


2.354 


















±0.068/ 


±0.016/ 


±0.036/ 


±0.016/ 














10 


-4.460 


2.277 


2.282 


2.313 


2.341 


2.352 
















±0.074/ 


±0.010/ 


±0.028/ 


±0.022/ 














11 


-4.498 


2.294 


2.287 


2.307 


2.331 


2.346 
















±0.074/ 


±0.010/ 


±0.019/ 


±0.022/ 


0.009/ 












12 


-4.518 


2.307 


2.293 


2.305 


2.324 


2.339 


2.345 














±0.072/ 


±0.017/ 


±0.011/ 


±0.019/ 


0.014/ 












13 


-4.540 


2.317 


2.299 


2.305 


2.319 


2.333 


2.341 














±0.068/ 


±0.022/ 


±0.004/ 


±0.015/ 


0.015/ 


±0.006/ 










14 


-4.559 


2.324 


2.305 


2.306 


2.316 


2.328 


2.337 


2.340 












±0.064/ 


±0.025/ 


±0.001/ 


±0.011/ 


0.014/ 


±0.009/ 










15 


-4.575 


2.330 


2.310 


2.308 


2.315 


2.324 


2.333 


2.338 












±0.061/ 


±0.027/ 


±0.005/ 


±0.007/ 


0.012/ 


±0.010/ 


±0.004/ 








16 


-4.589 


2.334 


2.315 


2.310 


2.314 


2.322 


2.330 


2.335 


2.338 










±0.057/ 


±0.027/ 


±0.008/ 


±0.004/ 


0.010/ 


±0.010/ 


±0.006/ 








17 


-4.602 


2.338 


2.318 


2.313 


2.314 


2.320 


2.327 


2.333 


2.336 










±0.053/ 


±0.027/ 


±0.010/ 


±0.001/ 


0.008/ 


±0.009/ 


±0.007/ 


±0.003/ 








-4.Dio 


o Q /I n 


Z.6ZI 


z.oio 


z.oio 


z.oiy 


z.ozo 


z.ooU 


O QQ /I 

z.oo4 


Z.OOD 








±0.050/ 


±0.027/ 


±0.012/ 


±0.001/ 


0.005/ 


±0.008/ 


±0.007/ 


±0.004/ 






19 


-4.623 


2.342 


2.324 


2.317 


2.316 


2.319 


2.324 


2.329 


2.332 


2.335 








±0.047/ 


±0.026/ 


±0.013/ 


±0.003/ 


0.003/ 


±0.007/ 


±0.007/ 


±0.005/ 


±0.002/ 




20 


-4.631 


2.344 


2.327 


2.319 


2.317 


2.319 


2.323 


2.327 


2.331 


2.333 


2.334 






±0.045/ 


±0.025/ 


±0.014/ 


±0.005/ 


0.002/ 


±0.005/ 


±0.006/ 


±0.007/ 


±0.003/ 





TABLE V: Comparison of Infrared and exact LDE at large- Ai', for r) = R^^ — l/(247r) at different orders k. 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




IR 


1.16424 


3.20165 


1.89188 


2.67411 


2.17385 


2.47173 


2.27958 


2.39236 


2.31806 


2.36035 




exact 


1.16424 


2.2129 


1.78534 


1.98289 


1.892 


1.94579 


1.92615 


1.94536 


1.94367 


1.95312 




k 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


100 


IR 


2.33129 


2.34685 


2.33524 


2.34073 


2.33590 


2.33765 


2.33549 


2.33588 


2.33481 


2.33474 


2.32910 


exact 


1.95614 


1.9624 


1.96651 


1.9715 


1.97566 


1.97999 


1.98393 


1.98782 


1.99149 


1.99504 


2.10987 
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TABLE VI; Best PMS results for ci, at different orders k, from the exact LDE series at large- A'^ when omitting the tadpole 
term -(A''77)/(47r). 



k 


3 


4 


5 


6 


10 


15 


Cl 


1.061 


1.222 ±0.37/ 


1.34 


1.435 ± 0.720/ 


1.710 ±0.99/ 


1.896 ± 1.089/ 


k 


20 


30 


50 


60 


80 


100 


Cl 


1.999 ± 1.120/ 


2.129 ±0.880/ 


2.243 ± 0.900/ 


2.272 ±0.910/ 


2.311 ±0.780/ 


2.33 ± 0.69/ 



TABLE VII: Standard LDE at large- A, Eq. limji . All CIRT-PMS resuhs for ci, at different orders k. 



k 


FO 


Fl 


F2 


F3 


F4 


F5 


F6 


F7 


F8 


F9 


FIO 


n 
Z 


-z.olo 


SI s 
z.olo 




















Q 

3 


O A An 

-3.4U0 


O AAA 

2.000 
























1 n Qi /I r 
±U.ol4i 




















4 


-3.004 


1 01 A 

1.219 


o /loo 
2.420 






















1 A C^iA 7" 




















r 




-3.09r 


A ^C^A 
0. <00 


O 1 oi 

2.121 










































6 


O 1 A 

-3.610 


A A£? 

0.006 


1 70 /I 

1. (34 


O O ,1 

2.240 




















1 A 0/1 7 r 


1 A Aao T 


















1 


-3.613 


A /I 1 

0.41D 


1 /I oo 

1.422 


O AAiJ 

2.1196 




















1 A 7/1 A r 


1 A QOfi r 


1 A OQQ r 
















o 
8 


-3.614 


A /I /I O 

U.448 


1 OAC; 

1.2U0 


1 ooo 
1.882 


O 1 /I 7 


















ztl. ZZOi 


1 A 1 ov r 

±U. 12 / 1 


1 A QA/1 r 
















9 


O tit A 

-3.614 


A K 07 

0.08/ 


1 A71 
1.0/ 1 


1 O 1 

1.6el 


O Af?0 

2.U62 


















11 f?Ai r 
ztzl.Dyli 


1 A 1 1 c r 
ztzU.llOJ 


1 A oflo r 


1 A A ACT 

rtU. 140-/ 














10 


-3.614 


0.836 


1.001 


1.518 


1.932 


2.087 
















1 o 1 QQ r 


1 A Qcc r 


1 A 1 71 r 


1 A nc\(i. T 
±U.2UDi 














11 


-3.614 


1 OAO 

1.203 


A AAO 

0.998 


1 O A7 

1.39 i 


1 7A7 

1. lyi 


O AO /I 

2.U34 
















Zt2.002i 


1 A f\f\f\ T 


1 A AQQ r 


1 A onn r 


1 A AA7 r 
zhU.Uy / 1 












1 o 


O HA A 

-3.614 


1 7A1 
1. (01 


1 A /I A 

1.04U 


1 1 7 

1.31 < 


1 /570 

1.6(8 


A A ,1 A 

1.949 


O A /I 7 

2.U4/ 














1 o QI Q r 

±2.y loi 


1 A Q/1 1 r 


1 A 1 1 1 r 
ztU. 1 1 IJ 


1 A 1 Cl r 
ztzU.lOlJ 


1 A AAOT 

±U. 14oJ 












13 


-3.614 


2.343 


1.126 


1.272 


1.581 


1.856 


2.012 














±3.190/ 


±1.054/ 


±0.266/ 


±0.073/ 


±0.148/ 


±0.067/ 










14 


-3.614 


3.142 


1.251 


1.58 


1.506 


1.768 


1.954 


2.020 












±3.339/ 


±1.243/ 


±0.419/ 


±0.021/ 


±0.122/ 


±0.102/ 










15 


-3.614 


4.106 


1.412 


1.272 


1.455 


1.691 


1.887 


1.996 












±3.301/ 


±1.401/ 


±0.566/ 


±0.125/ 


±0.073/ 


±0.110/ 


±0.048 








16 


-3.614 


5.229 


1.606 


1.309 


1.424 


1.628 


1.821 


1.955 


2.001 










±2.996/ 


±1.525/ 


±0.702/ 


±0.232/ 


±0.010/ 


±0.095/ 


±0.074 








17 


-3.614 


6.491 


1.828 


1.367 


1.411 


1.579 


1.761 


1.906 


1.984 










±2.327/ 


±1.604/ 


±0.825/ 


±0.339/ 


±0.062/ 


±0.063/ 


±0.082 


±0.035/ 






18 


-3.614 


7.838 


2.071 


1.442 


1.416 


1.544 


1.709 


1.855 


1.954 


1.988 








±1.169/ 


±1.632/ 


±0.932/ 


±0.442/ 


±0.140/ 


±0.019/ 


±0.073 


±0.055/ 






19 


-3.614 


9.181 


2.327 


1.532 


1.434 


1.522 


1.665 


1.807 


1.917 


1.975 








±0.628/ 


±1.600/ 


±1.022/ 


±0.539/ 


±0.220/ 


±0.033/ 


±0.051 


±0.061/ 


±0.025/ 




20 


-3.614 


10.360 


2.584 


1.634 


1.466 


1.513 


1.632 


1.764 


1.877 


1.952 


1.978 






±3.236/ 


±1.499/ 


±1.091/ 


±0.628/ 


±0.299/ 


±0.091/ 


±0.019 


±0.055/ 


±0.041/ 
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TABLE VIII: Infrared LDE at large- iV, Eq. KTEl . All CIRT-PMS large- iV results for ci, at different orders k. 



h 

K 


r i 


r L 


r o 


-1 

r 4 


r 


r D 


r I 


r o 


r y 


2 


2.328 


















3 


2.328 


2.262 


2.395 














4 


2.328 


2.320 


2.337 


— 






— 










±0.067/ 


±0.067/ 














5 


2.328 


2.369 


2.287 


2.271 


2.386 


— 


— 


— 


— 






±0.054/ 


±0.054/ 














6 


2.328 


2.389 


2.268 


2.294 


2.363 














±0.027/ 


±0.027/ 


±0.041/ 


±0.041/ 










7 


2.328 


2.391 


2.266 


2.323 


2.334 


2.281 


2.376 










±0.005/ 


±0.005/ 


±0.050/ 


±0.050/ 










Q 
O 


Z.oZo 


Z.OOD 


O 0^70 
1. 1 1 U 


O Q /I /I 


1.616 


z.zyo 


O OR A 
Z.OD4 










±0.011/ 


0.011/ 


±0.045/ 


±0.045/ 


±0.025/ 


±0.025/ 






9 


2.328 


2.379 


2.278 


2.357 


2.300 


2.310 


2.347 


2.289 


2.368 






±0.021/ 


±0.021/ 


±0.036/ 


±0.036/ 


±0.036/ 


±0.036/ 






10 


2.328 


2.371 


2.286 


2.363 


2.294 


2.325 


2.332 


2.296 


2.360 






±0.027/ 


±0.027/ 


±0.025/ 


±0.025/ 


±0.038/ 


±0.038/ 


±0.017/ 


±0.017/ 



TABLE IX: PMS optimization of Pade approximants (based both on the original and CIRT improved series for two of the 
relevant families) for ci in the large A*' case with the exact series coefficients. 



k 


P[l,2] PMS 




P[l,2] CIRT 




P[0,3] PMS 




P[0,3] CIRT 




3 


2.623 


2.034 


2.01 ± 0.14 I 




2.328 ± 0.18 I 




2.328 ± 0.051 




4 


2.61 


2.05 


1.877 ± 0.34 I 


2.02 


2.328 ± 0.17 I 




2.328 ± 0.051 




5 


2.615 ± 0.10 I 


2.04 ± 0.1 I 


1.77 ± 0.61 I 


1.98 ± 0.06 I 


2.29 ± 0.2 I 


2.37 ± 0.19 I 


2.34 ± 0.051 


2.33 ± 0.05 I 



TABLE X: All ordinary PMS and CIRT-improved PMS optimization for ci in the finite N = 2 case. 



k 


FO PMS 


FO CIRT 


Fl PMS 


Fl CIRT 


F2 PMS 


F2 CIRT 


2 


-3.05916 


-3.98590 


3.05916 


3.98590 






3 


-4.47035 


- 5.97078 


2.44730 


3.10543 












±1.65256/ 


±3.09300/ 






4 


-5.30592 


-7.03900 


1.53443 


1.19134 


3.14286 


5.22847 








±2.29581/ 


±4.33683/ 
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TABLE XI: Same as previous table for ci in the finite N = 2 case, but with IR large N perturbative coefficient Ki for i > 3 



h 


1/ u Jr IVIO 




pi PMC! 


pi PTRT 


po PMC! 


po PTRT 


2 


-3.059 


-3.986 


3.059 


3.986 






3 


-4.470 


- 5.971 


2.447 


3.105 












±1.653/ 


±3.093/ 






4 


-5.306 


-7.039 


1.534 


1.194 


3.143 


5.106 








±2.296/ 


±4.337/ 








71 7 






1 1 7fi 


'\ 71 










±2.83/ 


±4.33/ 






6 


-5.97 


-7.05 


1.29 


1.179 


4.00 


5.09 








±3.13/ 


±4.32/ 






10 


-6.43 


-7.05 


1.219 


1.179 


4.49 


5.09 








±3.66/ 


±4.32/ 







TABLE XII: PMS oi)tiinizatioii of Padc ajiproxiiriaiits (based both on the original and CIRT improved scries for two of the 
relevant families) for tiuit(> ,\' = 2 case. 



k 


P[l,2] PMS 




P[l,2] CIRT 




P[0, 3] PMS 




P[0, 3] CIRT 




3 


2.62 


2.04 


2.01 ± 0.20 I 




2.328 ± 3.24 I 




1.37 


3.28 


4 


2.60 


2.06 


1.98 ± 0.41 I 


1.98 


2.328 ± 3.07 I 




0.97 ± 0.34 1 


3.15 


5 


2.63 ± 0.07 I 


2.02 ± 0.07 I 


1.97 ± 0.62 I 


1.94 ± 0.12 I 


1.51 ± 3.45 I 


3.15 ± 3.45 I 


1.11 ± 0.44 I 


3.17 


10 


2.()() = 0.01 1 


2.05 = 0.01 1 


1.82 = 1 


1.88 = 0.11 1 


1.20 = l.l).") 1 


ii.l.") = !.().-) 1 


1.15 = 0.11 1 


:i.27 



